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1. INTRODUCTION 

1.1. What are electronic correlations and where do they show up? 

The term correlations (which means "with relation", from Latin con + relatio) is 
used not only in physics, but also in many other fields. For example, in grammar the two 
words either ... or are called a "correlate". Indeed, in a grammatically correct sentence 
the word either always has to be followed by the word or. Obviously the two words are 
correlated. In mathematics and natural sciences the term correlation is used to express 
the fact that the average or expectation value of a product of quantities is usually not 
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equal to the product of the averages of the individual quantities: 

(AB)?(A)(B). (1) 

For example, the mass or charge density n(r) of a many-body system at position r is 
influenced, in general, by the density of the particles at other positions r' . Therefore the 
density-density correlation function is not simply given by the product of the average 
density n: 

(n(r)n(r>))^(n(r)}(n(r>)}=n 2 . (2) 

Indeed, this very property defines correlations: they express effects which go beyond 
approximations obtained by the factorization of correlation functions, i.e., beyond static 
mean-field theories such as the Weiss mean-field theory for the Ising model or the 
Hartree approximation for the Hubbard model. 

Correlations in space and time are by no means abstract notions, but occur frequently 
in everyday life. Persons in an elevator or in a car are strongly correlated both in space 
and time, and it would be quite inadequate to describe the situation of a person in such 
a case within a factorization approximation where the influence of the other person(s) is 
described only by a static mean-field, i.e., a structureless cloud. 

As in the case of two persons riding together on an elevator, two electrons with 
different spin direction occupying the same narrow d or / orbital in a real material 
are also correlated. Here the degree of correlation can be estimated in a very simplified 
picture as follows. Assuming the correlated electrons (or rather the quasiparticles, i.e., 
excitations) to have a well-defined dispersion % their velocity is given by = \ | V^e^l . 
The typical velocity is given by ~ f , where a is the lattice spacing and x is the 
average time spent on an atom. The derivative can be estimated as \ | Vk^kl ~ \dW since 
|Vjt| ~ l/k ~ a and |gjt| corresponds to the band overlap t and hence to the band width 
W. Altogether this means that 



The narrower an orbital, the longer an electron therefore resides on an atom and thereby 
feels the presence of other electrons. Hence a narrow band width implies strong elec- 
tronic correlations. 

Indeed this is the case for many elements in the periodic table. Namely, in many ma- 
terials with partially filled d and / electron shells, such as the transition metals V, Fe, 
and Ni and their oxides, or rare-earth metals such as Ce, electrons occupy narrow or- 
bitals. This spatial confinement enhances the effect of the Coulomb interaction between 
the electrons, making them "strongly correlated". Correlation effects lead to profound 
quantitative and qualitative changes of the physical properties of electronic systems as 
compared to non-interacting particles. In particular, they often respond very strongly 
to changes in external parameters. This is expressed by large renormalizations of the 
response functions of the system, e.g., of the spin susceptibility and the charge com- 
pressibility. Electronic correlations also play an essential role in high temperature super- 
conductivity. In particular, the interplay between the spin, charge and orbital degrees of 
freedom of the correlated d and / electrons and with the lattice degrees of freedom leads 
to a wealth of unusual phenomena at low temperatures [1]. These properties cannot be 
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explained within conventional mean-field theories, e.g., Hartree-Fock theory, since they 
describe the interaction only in an average way. 



1.2. Electrons vs. Landau quasiparticles 

Electrons are fermions which obey the Fermi-Dirac statistics. The Pauli exclusion 
principle implies the existence of a Fermi body of occupied states and thereby of a Fermi 
surface which distinguishes between occupied states (inside the Fermi body) and empty 
states (outside the Fermi body). The existence of a Fermi surface is something quite 
extraordinary. It allows for the formulation of Landau Fermi-liquid theory and thereby 
for a deep understanding of interacting fermionic systems [2]. 

In the Landau Fermi-liquid theory a one-to-one correspondence between fc-states of 
the non-interacting and the interacting system is assumed. Therefore there exist well- 
defined ^-states, called quasiparticles, which have a finite lifetime (the closer they are to 
the Fermi surface the more well-defined they are, i.e., the longer they live), an effective 
mass m* and an effective interaction. Quasiparticles are the elementary excitations of 
a Fermi liquid and determine the entire low temperature thermodynamics of a Landau 
Fermi liquid. They are rather abstract objects which should not be confused with the 
particles of the non-interacting system. The Fermi-liquid concept is independent of the 
strength of the bare interaction between particles. Hence it can not only describe simple 
metals such as potassium where m* is not much different from the bare electronic mass, 
but even "heavy fermion systems" such as UBei3 where m* can be a factor 1000 larger 
than the bare electronic mass [3]. 



1.3. The simplest model for correlated electrons 

The simplest model describing interacting electrons in a solid is the one-band, spin- 
1/2 Hubbard model [4-6] where the interaction between the electrons is assumed to be 
so strongly screened that it is taken as purely local. The Hamiltonian consists of two 
terms, the kinetic energy Ho and the interaction energy H[ (here and in the following 
operators are denoted by a hat): 

H = Ho + Hj (4a) 

#o = EE^ Jff = E^ (4b) 
Ri,Rj a k,a 

H l = ( 4c ) 

Ri 

where cf a (c ia ) are creation (annihilation) operators of electrons with spin o at site Ri, 
and h ia = c^ a c ia . The Fourier transform of the kinetic energy in (4b), where tjj is the 
hopping amplitude, involves the dispersion £ k and the momentum distribution operator 

"to- 
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FIGURE 1. Schematic illustration of interacting electrons in a solid in terms of the Hubbard model. The 
ions appear only as a rigid lattice (here represented as a square lattice). The electrons, which have a mass, 
a negative charge, and a spin (f or I), move from one lattice site to the next with a hopping amplitude 
t . The quantum dynamics thus leads to fluctuations in the occupation of the lattice sites as indicated by 
the time sequence. When two electrons meet on a lattice site (which is only possible if they have opposite 
spin because of the Pauli exclusion principle) they encounter an interaction U . A lattice site can either be 
unoccupied, singly occupied ( j" or J,), or doubly occupied. 



A schematic picture of the Hubbard model is shown in Fig. 1. When we look only at 
a single site of this lattice model, this site will sometimes be empty, singly occupied or 
doubly occupied. In particular, for strong repulsion U double occupations are energet- 
ically very unfavorable and are therefore strongly suppressed. In this situation (nj+n^) 
must not be factorized since (n^n^) ^ (n^)(nii). Otherwise, correlation phenomena 
such as the Mott-Hubbard metal-insulator transition are eliminated from the very begin- 
ning. This explains why Hartree-Fock-type mean-field theories are generally insufficient 
to explain the physics of electrons in the paramagnetic phase for strong interactions. 

The Hubbard model looks very simple. However, the competition between the kinetic 
energy and the interaction leads to a complicated many-body problem, which is impos- 
sible to solve analytically, except in dimension d = 1 [7]. This model provides the basis 
for most of the theoretical research on correlated electrons during the last decades. 

2. MEAN-FIELD THEORIES FOR MANY-BODY SYSTEMS 

2.1. Construction of mean-field theories 

It is well known that theoretical investigations of quantum-mechanical many-body 
systems are faced with severe technical problems, particularly in those dimensions 
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which are most interesting to us, i. e., d = 2, 3. This is due to the complicated dynamics 
and, in the case of fermions, the non-trivial algebra introduced by the Pauli exclusion 
principle. In the absence of exact methods there is clearly a great need for reliable, 
controlled approximation schemes. Their construction is not straightforward. 

In the statistical theory of classical and quantum-mechanical systems a rough, overall 
description of the properties of a model is often obtained within a so-called mean-field 
theory. Although the term is frequently used, the actual meaning of what a mean-field 
theory is or should be is rather vague, because there is no unique prescription of how to 
construct such a theory (and sometimes it is only a respectable name for a questionable 
approximation . . . ). Hence every time one encounters the term "mean-field theory", one 
should ask about the reliability of this approximation, i. e., about its range of validity 
with respect to the input parameters and its thermodynamic consistency. 

There exists a well-established branch of approximation techniques which makes 
use of the simplifications that occur when some parameter is taken to be large (in 
fact, infinite), e.g., the length of the spins 5, the spin degeneracy N, the coordination 
number Z (the number of nearest neighbors of a lattice site). Investigations in this limit, 
supplemented, if at all possible, by an expansion in the inverse of the large parameter, 
often provide valuable insight into the fundamental properties of a system even when 
this parameter is not large. 

One of the best-known mean-field theories is the Weiss molecular-field theory for the 
Ising model [8]. It is a prototypical single-site mean-field theory which becomes exact 
for infinite-range interaction, as well as in the limit of the coordination number Z — > <*> 
or 1 the dimension d — > °°. In the latter case 1/Zorl/Jisa small parameter which 
can sometimes be used to improve the mean-field theory systematically. This mean- 
field theory contains no unphysical singularities and is applicable for all values of the 
input parameters, i.e., coupling parameters, magnetic field, and temperature. It is also 
diagrammatically controlled [10]. Insofar it is a very respectable approximation which 
sets very high standards for other mean-field theories. 

2.7.7. Motivation for using the limit of high dimensions to construct 

mean-field theories 

In a perfectly crystalline system every lattice site has the same number of nearest 
neighbors Z. In three dimensions (d = 3) one has Z = 6 for a simple cubic lattice (Z = 2d 
for a hypercubic lattice in general dimensions d), Z = 8 for a bcc lattice and Z = 12 for an 
fcc-lattice. The dimensionality of a lattice system is directly described by the number Z 



For regular lattices, e.g., Bravais-lattices, both a dimension d and a coordination number Z can be 
defined. In this case either d or Z can be used alternatively as an expansion parameter. However, there 
exist other lattices (or rather graphs) which cannot be associated with a physical dimension d although a 
coordination number Z is well-defined. The best-known example is the Bethe lattice, an infinitely extended 
Cayley tree [8, 9], which is not a regular lattice because it does not have loops. The coordination number Z 
is therefore a very useful parameter for theoretical investigations of lattice models, although the dimension 
d is the more general physical parameter. In the following discussion we mostly use both d and Z in 
parallel. 
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(see footnote 1). Since Z ~ <^(10) is already quite large in d = 3, such that l/Z is rather 
small, it is only natural and in the general spirit of theoretical physics to consider the 
extreme limit Z — > oo to simplify the problem. Later, if possible, one can try to improve 
the result by expanding in the small parameter l/Z. The limit d — > °° is not as academic 
as it might seem. In fact, it turns out that several standard approximation schemes which 
are commonly used to explain experimental results in dimension d = 3, are exact only 
in d = oo [11]. 



2.2. A prototypical example: The Weiss mean-field theory for the 

Ising model 

In the case of classical spin models (e.g., Ising, Heisenberg) the Z — > oo limit is well- 
known [8, 10]. It leads to the results of the Weiss molecular-field theory which may be 
viewed as the prototypical method for constructing a mean-field theory. The Hamiltonian 
for the Ising model with nearest-neighbor (NN) coupling is given by 

H=-\j £ SiSj, (5) 

where we assume ferromagnetic coupling (J > 0). Every spin 5,- interacts with a local 
field hi, produced by its nearest neighbors at site In the Weiss mean-field approach the 
two-spin interaction in (5) is decoupled, i. e., H is replaced by a mean-field Hamiltonian 

H MF = -h MF J £S l +E shm . (6a) 
Ri 

Now a spin Si interacts only with a global ("molecular") field 

(0 

h M F = /£(5y> (6b) 

= J(S). (6c) 

Here ( ) indicates the thermal average, Eshift = \ LJZ(S} 2 is a constant energy shift with 
L as the number of lattice sites and the superscript (?) implies summation over only 
NN-sites of Ri. This corresponds to the factorization 

([ 5i --(5)][5 y -(5)]> = 0, (7) 

whereby correlated fluctuations of spins at sites Ri and Rj are neglected. In the limit 
Z — > oo the coupling constant J has to be rescaled as 

J* 

J ->■ — , J* = const (8) 

z 

for /zmf to remain finite. In this limit the factorization procedure (7), and hence the 
replacement of (5), by the mean-field Hamiltonian (6a), becomes exact [12, 13]. 
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FIGURE 2. Already in three dimensions (d = 3) the coordination number Z of a lattice can be quite 
high, as in the face-centered cubic lattice where Z = 12. In the limit Z — > oo, or equivalently d — > °°, the 
Ising model effectively reduces to a single-site problem where the local field /i, is replaced by a global 
mean ("molecular") field h^e- 



Eq. (6a) implies that in the limit Z — > oo fluctuations in the "bath" of surrounding 
neighbors become unimportant, such that the surrounding of any site is completely 
described by a single mean-field parameter /?mf (see Fig. 2). Hence the Hamiltonian 
becomes purely local 

H MF = ^H l+ E shm (9) 
Ri 

Hi = -hMpSi. (10) 

Thereby the problem reduces to an effective single-site problem. The value of (S) is 
determined by the Curie-Weiss (or Bragg-Williams) self-consistent equation 

(S)=tanh(j8J*(5)), (11) 

where /3 = 1 /T (here kg = 1 ) . 

It should be noted that the scaling (8) is typical for localized spin models with 
isotropic coupling, i. e., when the spatial average = J is non-zero. On the other hand, 

in the classical spin glass problem with random coupling one has Jij = 0, but Jfj ^ 0. In 

this case a different kind of scaling, namely Jjj — >■ JfJ yZ, has to be used [14]. 

2.3. Hartree approximation for the Hubbard model 

To elucidate some of the shortcomings of conventional mean-field theories for 
Hubbard-type models, we consider the Hartree approximation 2 , which corresponds to a 
factorization of the interaction term in direct analogy to the factorization of the spins in 



The Hubbard interaction acts only between electrons with opposite spin; therefore the Hartree-Fock 
approximation does not lead to an exchange term. 
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the Ising model. To this end let us consider a generalization of Hubbard model 

H = H^+\ £ £ Vifn^, (12) 
Z Ri,RjOo' 

where the interaction V^ CT ' describes a local part, U8 a - a i, i. e., the Hubbard interaction, 

as well as a nearest-neighbor contribution V aa ' . The form of the interaction term is sim- 
ilar to (5) for the Ising model. In the spirit of a conventional mean-field approximation 
the two-particle interaction in (12) is factorized, i. e., H ext is replaced by 

# MF = #kin + £ fi ia (h ia ) +£shift (13) 

in complete analogy with the Weiss molecular field theory for the Ising model in 
Sec. 2.2. In (13) a a-electron at site /?, interacts only with a local field (a c-number) 

(0 

(hie) = £ V? a (h ja ). (14) 

The above decoupling of the operators is equivalent to the Hartree approximation which 
sets 

([bio- (n ia )][njo' -{*]&)]) =0 (15) 

for all a, <r'. Thereby correlated fluctuations on the sites /? ; and /? 7 are neglected. 
Although (13) is a one-particle problem it cannot be solved exactly in any systematic 
way, since, in principle, the potential, i. e., the mean field (hi a ), may be an arbitrarily 
complicated function of position. 

An obvious question concerning (13) is, whether the mean-field decoupling ever 
becomes exact for all values of U,V aa and n, i. e., beyond the weak-coupling or 
low-density limit. How about the limit d — > °°? We will come back to this question 
in Sec. 3.3.1 once we understood how this limit can be employed in the case of lattice 
fermion models. 

3. LATTICE FERMIONS IN THE LIMIT OF HIGH DIMENSIONS 

3.1. Scaling of the hopping amplitude 

It is natural to ask whether the limit d — > °° may also be useful in the investigation of 
lattice models with itinerant quantum-mechanical degrees of freedom and, in particular, 
in the case of the Hubbard model. Following Metzner and Vollhardt [15] we take a look 
at the kinetic energy term (4b), since the interaction term is purely local and is thereby 
completely independent of the lattice structure and the dimension. For nearest-neighbor 
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hopping on a d-dimensional hypercubic lattice (where Z = 2d) is given by 3 



£fc = —2t^ coski. 
i=l 

The density of states (DOS) corresponding to is 

N d (co) =£5(^w -e k ). 
k 



(16) 



(17) 



This is simply the probability density for finding CO = for a random choice of k — 
(ki,.. .,kd). If the ki are chosen randomly, in (16) is the sum of (independent) random 
numbers —2t coski. The central limit theorem then implies that in the limit d — > °° the 
DOS is given by a Gaussian 



N d (co) 



1 



2tV7td 



exp 



(-Y 



(18) 



Unless t is scaled properly with d this DOS will become arbitrarily broad and featureless 
for d — > °o. Clearly only the "quantum" scaling 



t — > —=, t = const., 
yd 



(19) 



yields a non-trivial DOS [16], [15]: 

Noo{(6) = 



2m* 



exp 



1 /co\ 2 

2 W 



(20) 



This DOS does not have any van Hove singularities. The reason for this can be seen 
when Nd(co) is calculated explicitly from (17) [16, 17]. Expressing the 5-function as a 
Fourier series one has 



h{co-e k ) 



(21a) 



= / Jre 



(COT [ 



/o(2«) 



(21b) 



where 7oW = 1 — * 2 + @{x A ), x <ti 1, is the zero-order Bessel-function. For d 3> 1 the 
main contribution to the integral comes from the first extremum of Jq(x), i. e., | x |~ 



3 In the following we set Planck's constant Ti, Boltzmann's constant kg, and the lattice spacing equal to 
unity. 
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FIGURE 3. Density of states of a tight-binding dispersion relation on a hypercubic lattice in d 
1,2,3,4,5 (full lines) as compared to the result for d = °° (dashed line); from Ref. [11]. 



\/2ty/d, while van Hove singularities are due to higher extrema, yielding exponentially 
small contributions to N c i((o). Hence, using the scaling (19), one finds for d 3> 1 



N d {(o) = 



2k t* 



h(o>in 2 \ !__L 



16 J 



(f) 4 -(f) 2 + 3 



(?) 



(22) 



It is interesting to compare Nd(co) for different d as shown in Fig. 3. For d > 3 the shapes 
rapidly approach the d = °° result; the main difference is that for d < °° the band has a 
finite width, while in d = °° there exist exponentially small tails for all ft). 
The scaling (19) expresses the fact that, for some randomly chosen k, 



e k /t ~ &{Vd) , d 



(23) 
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since e^/t is the sum of d — > °° many random numbers from the interval [—1,1]. Clearly, 
k = and k = {k, . . . , 71) are special values, for which (23) does not hold. However, as 
long as £fc appears under an integral these points have zero measure. 
Using the relation 

t£^(£*)= / dcoN(co)F(co) (24) 

k J —<x> 

where L is the number of lattice sites, the kinetic energy of the non-interacting electrons 
is found as 

EL = £ Who = -2Lf 2 N„{E F ). (25) 

We see that only the quantum scaling (19) leads to a finite kinetic energy in d = °°. 
Without this scaling E^ n would diverge, while the classical scaling t —¥ t*/Z would 
reduce to zero 4 . 

The interaction term in (4) is seen to be purely local and independent of the surround- 
ing; hence it is independent of the spatial dimension of the system. Consequently, the 
on-site interaction U need not be scaled. So we see that the scaled Hubbard Hamiltonian 

7f E Ef+c> + (/£« lt % (26) 
(Ri.Rj) o Ri 

has a nontrivial Z — > <*> limit, where both terms, the kinetic energy and the interaction, 
are of the same order of magnitude and are thereby able to compete. It is this competition 
between the two terms which leads to interesting many-body physics (see footnote 4). 

The quantum scaling (19) was determined within a k- space formulation. We will now 
derive the same result within a position-space formulation. 



3.2. Simplifications of the many-body perturbation theory 

The most important consequence of the scaling (19) is the fact that it leads to signifi- 
cant simplifications in the investigation of Hubbard-type lattice models [15, 17, 19-22]. 



4 To obtain a physically meaningful mean-field theory for a model its internal or free energy has to remain 
finite in the limit dorZ^oo. While for the Ising model the scaling J — > J/Z, J= const., was rather 
obvious this is not so for more complicated models. Namely, fermionic or bosonic many-particle systems 
are usually described by a Hamiltonian consisting of several non-commuting terms, e.g., a kinetic energy 
and an interaction, each of which is associated with a coupling parameter, usually a hopping amplitude 
and an interaction, respectively. In such a case the question of how to scale these parameters has no unique 
answer since this depends on the physical effects one wishes to explore. In any case, the scaling should 
be performed such that the model remains non-trivial and that its internal or free energy stays finite in the 
Z — > °o limit. By "non-trivial" we mean that not only (Hq) and (Hi nt ), but also the competition between 
these terms, expressed by ([Ho, Ant]), should remain finite. In the case of the Hubbard model it would be 
possible to employ classical scaling for the hopping amplitude, i.e., t — > t* /Z, t* = const., but then the 
kinetic energy would be reduced to zero in the limit d — > °°, making the resulting model uninteresting (but 
not unphysical) for most purposes. For the bosonic Hubbard model the situation is more subtle due to the 
occurrence of Bose-Einstein condensation; for a discussion see Ref. [18]. 
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To understand this point better we take a look at the perturbation theory in terms of U. 
At T = and U = the kinetic energy (25) may be written as 

EL=~t £ £ Sj.cn (27) 

where g9. CT = {cf a 6j a )o is the one-particle density matrix. This quantity can also be 
interpreted as the amplitude for transitions between site /?, and Rj, whose square is 
proportional to the probability for a particle to hop from Rj to i.e., | #9. ^ | 2 ~ 1/Z ~ 
l/J since /?j has <^(d) nearest neighbors J?,-. Thus the sum of | gfj a | 2 over all nearest 
neighbors must yield a constant. In the limit rf->»we then have 

4. G ~^(-^),tf ? AWoftf;. (28) 

Since the sum over the NN-sites /?, in (27) is of ^(d) the NN-hopping amplitude t must 
obviously be scaled according to (19) for E^ n to remain finite in the limit J,Z — > °«. 
Hence, as expected, a real-space formulation yields the same results for the required 
scaling of the hopping amplitude. 

The one-particle Green function ("propagator") G?. a ((Q) of the non-interacting sys- 
tem obeys the same scaling as CT . This follows directly from its definition 

G^ a (t) = -(Tc ia (t)c+ a (0)) , (29) 

where T is the time ordering operator, and the time evolution of the operators is provided 
by the Heisenberg representation. The one-particle density matrix is obtained as gfj a = 

lim^o- G®j a (t). If g®j a obeys (28) the one-particle Green function must follow the 
same scaling at all times since this property does not dependent on the time evolution and 
the quantum mechanical representation. The Fourier transform G^. a {(6) also preserves 
this property. 

It is important to realize that, although the propagator G°- a ~ 1 / \fd vanishes for 
d — > oo, the particles are not localized, but are still mobile. Indeed, even in the limit 
d — > °° the off-diagonal elements of G°. a contribute, since a particle may hop to d 

nearest neighbors with reduced amplitude t*/y/2d. For general i, j one finds [19, 23] 

Gfj^^tffl/JK-*^, (30) 

where || if ||= Lf=i I Rn I is the length of R in the so-called "New York metric" (also 
called "taxi cab metric", since particles only hop along horizontal or vertical lines, never 
along a diagonal). 

It is the property (30) which is the origin of all simplifications arising in the limit 
d — > oo. In particular, it implies the collapse of all connected, irreducible perturbation 
theory diagrams in position space [15, 17, 19]. This is illustrated in Fig. 4, where a 

(2) 

contribution in second-order perturbation theory to the irreducible self-energy, E - • , is 
shown. 
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FIGURE 4. Contribution to the irreducible self-energy for the Hubbard model in second-order pertur- 
bation theory in U, and its collapse in the limit d -> °°. 



In all physically interesting circumstances E>y will only enter in a sum over /?, and Rj. 
Then it becomes apparent that E^ is only of order 1 / \fd small, unless i = j. Namely, 

for j 7^ i the three lines, corresponding to G°- a , contribute a factor 1/d 3 / 2 , while the 
sum over NN- sites Rj of /?, contributes a factor J. Only for i = j is the value of 

(2) 

L]j J independent of d. Hence in the limit d — > °° the diagram on the left-hand side of 
Fig. 4 is equivalent to the "collapsed", petal-shaped diagram on the r.h.s., provided i = j; 
otherwise it is zero. More generally, any two vertices which are connected by more than 
two separate paths will collapse onto the same site 5 . In particular, the external vertices 
of any irreducible self-energy diagram are always connected by three separate paths and 
hence always collapse. As a consequence the full, irreducible self-energy becomes a 
purely local quantity [15, 17]: 

E \a(«) rf =~£«,a(<»)S i7 . (31a) 

In the paramagnetic phase we may write E, !i(T ((o) =£(©). The Fourier transform of 
Ejj fj is seen to become momentum-independent 

E f7 (*,(0) d = oo E f7 ((0). (31b) 

This leads to tremendous simplifications in all many-body calculations for the Hubbard 
model and related models. It should be noted that a ^-independence of E is sometimes 
assumed as a convenient approximation ("local approximation") [24-26]. Here we iden- 
tified the limit where this is indeed exact. 

The result expressed in (3 lb) may equally be obtained by working in &-space from 
the beginning [17]. For this we consider an external vertex where a momentum k enters 
from outside (see Fig. 5). 

Making use of the fact that 
(i) the Hubbard interaction is momentum independent, 



5 Here a "path" is any sequence of lines in a diagram; they are "separate" when they have no lines in 
common. 
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k :J = k + k 2 - ki 



FIGURE 5. Interaction vertex for the Hubbard model with external momentum and frequency k 
(*,©). 



(ii) the ^-dependent, free propagator is given by 



1 



(d-£ k + n 



G°(e*,<») 



(32) 



and only depends on k via £ k (we do not write the spin index explicitly), and 

(iii) that momenta only enter explicitly in the conservation of momentum at the vertex, 

we see that the evaluation of the vertex involves the momentum summation 



1_ £ Gl((D 1 )Gl((K)Gl(c0 3 )8*(k-k 1 +k2-k 3 ) 

3 r r i 

Yl / d£iG°(£i, oh) N k (ei , e 2 , £ 3 ) , 
i=i lJ J 



where C0\ — Qh_ + Oh, = 00 and 



Jq-(R-R Q ) 



(33a) 
(33b) 

(34) 



is the "Laue-function" which guarantees momentum conservation up to a reciprocal 
lattice vector K. The lattice summation extends over all sites R relative to some origin 
Rq (without loss of generality we may put Rq = 0) and couples momenta explicitly. 
Without this coupling the momentum- summation would be simple because we would 
be able to use (24). In (33b) we therefore introduced a generalized density of states 

N k (£ h £ 2 ,£ 3 ) = ^ £ S(£ l -£ kl )8(£ 2 -£ k2 )8(£3-e k3 )S*(k-k l +k 2 -k 3 ), (35) 

k x ,k 2 M 

which is the probability density for £ ; = £^.(i = 1,2,3) for given k. Writing the <5- 
functions in (35) as a Fourier series (see (21a)) and using (34), one finds in the limit 

d = 00 

N k (£ h £2,£3)=N 00 (£i)N 00 (£ 2 )Noo(£3), d = 00 (36) 

i. e., for all k the quantity N k factorizes into a product of one-particle DOS's. This is 
equivalent to replacing the Laue-function in (35) by unity 



8*(q) = 1. 



(37) 
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It effectively means that the momentum conservation constraint may be ignored in 
d = oo. Defining the position of the interaction vertex in Fig. 5 by Rq = 0, (37) means 
that in the lattice sum over R only the local term R = contributes. This is, once 
again, the collapse-phenomenon discussed above. Due to the irrelevance of momentum 
conservation an external momentum k cannot enter into the internal structure of an 
irreducible self-energy diagram; this makes the irreducible self-energy ^-independent 
in d = oo (see (31b)). Note, however, that the total momentum of a particle must be 
conserved for the theory to be meaningful. 

Due to the simplifications caused by (36) or (31), the most important obstacle for 
actual diagrammatic calculations in finite dimensions d > 1, namely the integration 
over intermediate momenta, is removed in d = °°. While in finite dimensions these 
integrations lead to untractable technical problems, they become simple in d = oo, since 
one can replace them by one-dimensional integrations over the DOS. 

It should be noted that the limit d — > oo does not affect the dynamics of the system at 
all. Time is always one-dimensional and hence there is no "collapse" in the frequency 
variables. In spite of the simplifications in position (or momentum) space the problem 
retains its full dynamics in d = oo. 



3.3. Interactions beyond the on-site interaction 

In the case of more general interactions than the Hubbard interaction, e. g., nearest 
neighbor interactions such as 

Ann = £ £ V aa m ia n ja , (38) 
(Ri,Rj) era' 

the interaction constant has to be scaled, too, in the limit d — > oo. In the case of (38), 
which has the form of a classical interaction, the "classical" scaling 

V* , 

V aa > -+ -f- (39) 

is required. Of course, the propagator still has the dependence (30). The self-energy has 
the general diagrammatic form shown in Fig. 6. 

Due to (39) all contributions, except for the Hartree-term, are found to vanish in 
d = oo. Hence nonlocal interactions only contribute via their Hartree-contribution, which 
is purely static. This gives the Hubbard interaction a unique role: of all interactions for 
fermionic lattice models only the Hubbard interaction remains dynamical in the limit 
d-X»[17]. 



3.3.1. The Hartree approximation revisited 

We are now able to answer the question whether the decoupling leading to the Hartree 
approximation (13) also becomes exact in the limit d — > oo as in the case of the Ising 
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0(1) 0(l/v/d) 0(l/Vd) 

FIGURE 6. Diagrams representing the irreducible self-energy. 



model. Using (15) the answer is simple: as long as the on-site interaction in (12) is 
present, the decoupling (15) cannot become exact in d — > °°, because for i = j (and 
hence o = —o') the r.h.s. of (15) is always of order unity. Namely, at a given site Ri 
the potential, in units of U, felt by an electron with spin a is either (if there is no 
— o spin present) or 1 (if a —a spin is present). These on-site fluctuation effects are not 
described by the Hartree decoupling approximation. This has been explicitly verified by 
investigating the thermodynamics of (12) in d = °° [27]. On the other hand, if the on- 
site interaction in (12) did not exist (e.g., in the case of spinless fermions) the Hartree 
approximation (15) would indeed become exact in the limit d — > oo. 



3.4. One-particle and two-particle propagators 

Due to the ^-independence of the irreducible self-energy, (31b), the one-particle 
propagator of an interacting lattice fermion system is given by 

G ka (G)) = — jr-r^r- ( 4 °) 

Most importantly, the k dependence of G^co) comes entirely from the energy dispersion 
££ of the non-interacting particles. This means that for a homogeneous system with the 
propagator 

G ihc {(D)=L- l Y,GkA G >y HRi ~ Ri) ( 41 ) 
k 
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FIGURE 7. Diagrammatic contribution to the two-particle propagator. 



its local part, i. e., Gu, a , can be calculated in closed form [20] 

oo 

Woo(E) 



G i/)ff (fi>) = L- 1 £G )kia ((0)= / J£ 



(0-£' + /i-E f7 ((o) 



(42a) 



|_ e 4 e rfc(-/z CJ ) 



(42b) 



(42c) 



where z a = (co + /i - E CT (cu))/\/2?*, N^E) is given by (20), and erfc(x) is the error 
function. In the following we will limit our discussion to the paramagnetic phase and 
omit the spin index. The spectral function of the interacting system (often referred to as 
the DOS as in the non-interacting case) is then given by 



A(ffl) = --ImG(o) + /0+); 
it 



(42d) 



for U = one has A (to) = N((Q). In the limit d — > °° two quantities then play the most 
important role: the local propagator G{(6) and the self-energy E(co). 

Concerning the two-particle propagator G2 (or correlation functions, etc.) the collapse 
phenomenon is a little different from the one discussed below (35). Namely, it does not 
occur for all external momenta [23, 28]. A typical contribution to G2 is shown in Fig. 7, 
where Ti,r2 are irreducible, momentum independent vertices. 

The calculation of this contribution is analogous to that in (33), and involves the k- 
sum over a product of two one-particle propagators, i. e., over with k 2 = k\ + k. 

Introducing a generalized DOS 



N k (e h e 2 ) = £ 8(£i-£ kl )8(£2-e k2 )8*(k 2 -k 1 -k) 



(43) 



we find that for d = °° this quantity only factorizes into A^(ei) Ncofa) if k 7^ 0, 
k^ (k,...,k), i. e., when £ k ~ &(l/y/d). Hence the Bethe-Salpeter ladder collapses 
for all k, except for these two special values. 
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FIGURE 8. Correlation energy for the Hubbard model in second-order perturbation theory in U, 

e-i = j (2U 2 / | Co | ), vs. density n for several lattice dimensions. Left: d = 1,2,°°; right: d = 3,°°. 
Here | £ | is the kinetic energy for U — and n = 1; after Ref. [15]. 

3.5. Weak-coupling correlation energy for the Hubbard Model 

The correlation energy E C (U) is defined as the energy by which the Hartree-Fock 
energy is lowered when genuine correlations are included 

E C (U) := E exact (U)-E HF (U). (44) 

For d — > oo second-order Goldstone perturbation theory in the weak interaction U yields 
the simple expression (n^ = n^ = \) [15] 

eP(U) = -U 2 [°° d(oe 2c ° 2 P 2 (E F -(o)P 2 (-E F -(o), (45) 
Jo 

where P(x) is the Gaussian probability function. It should be stressed that for general 

(2) 

dimensions the calculation of E c involves 3d momentum integrals over a singular 
integrand. Analytic calculations are thereby ruled out and even numerical integration 
techniques becomes very cumbersome. Indeed, the application of Monte-Carlo inte- 
gration techniques becomes mandatory already for d ~ 2. By contrast, the case d = °° 
is seen to be the simplest of all dimensions, including d = 1, since it only demands a 
one-dimensional integral. 

In Fig. 8 E^ j V 2 is shown as a function of particle density n for several dimensions 
d. We see that the result for d = 3, which can only be obtained by considerable effort, 
is very well approximated by the result for d = °°, which is easily calculated. The 
differences between d = 1 , 2 and d = 3 are mainly due to the singularities of the DOS in 
d = l,2. 
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3.6. Consequences of the ^-independence of the self-energy 



We now discuss some more consequences of the ^-independence of the self-energy as 
derived by Muller-Hartmann [20]. Let us consider the Hubbard model and concentrate 
on the paramagnetic phase. That is, we do not wish to discuss the more complicated 
situation with a broken symmetry for the moment. 6 In the paramagnetic case at T = 
the one-particle propagator (40) takes the form (again we do not write the spin index 
explicitly) 

Gk(co) = =7-^-. (46) 

In general, even when E(ft>) is ^-dependent, the Fermi surface is defined by the ft) = 
limit of the denominator of (46) as 

£ k + L k (0)=E F . (47a) 

According to Luttinger and Ward [29] the volume within the Fermi surface is not 
changed by interactions, provided the effect of the latter can be treated in infinite-order 
perturbation theory (hence no broken symmetry). This is expressed by 

n = Y,e[E F -£ k -L k (0)}, (47b) 

where n is the particle density and 6(x) is the step function. In general, the ^-dependence 
of £fc(0) in (47a) implies that, in spite of (47b), the shape of the Fermi surface of the 
interacting system will be quite different from that of the non-interacting system (except 
for the rotationally invariant case £% = /(|Jfc|). For lattice fermion models in d = °°, 
where Sjt(co) = £(ft)), (47a) implies that the Fermi surface itself (and hence the volume 
enclosed) is not changed by interactions. 7 The Fermi energy is simply shifted uniformly 
from its non-interacting value E F , i. e., Ep = E F + E(0), to keep n in (47b) constant. 
From (42b) we thus conclude that the ft) = value of the local propagator, G(0), and 
hence of the spectral function, A(0) = — ^ImG(/0 + ), is not changed by interactions. 8 
Renormalizations of N(0) can only come from a ^-dependence of E, i. e., if dZ/dk^ 0. 
For ft) — > the self-energy has the property 

ImE(ft>)°c(o 2 (47c) 

which implies quasiparticle (Fermi liquid) behavior. The effective mass 



m* _ i dL 
m dCO 



1 r , Im£(ft) + z0~) , „ , 

= 1 + - / dCO > 1 (47d) 

K J-oo CO 2 



6 On bipartite lattices and for nearest neighbor hopping the Hubbard model has a "perfect-nesting insta- 
bility" at half filling due to ^ k = —E^+q with Q = (n,...,n), leading to an insulating state with antiferro- 
magnetic correlations in d > 2. In this case the system is not in a paramagnetic state even for arbitrarily 
small U. However, by including a finite hopping amplitude to next-nearest neighbors the paramagnetic 
phase becomes stable for small U. 

1 In d — °° the notion of a Fermi surface of a lattice system is made complicated by the fact that in this 
limit the dispersion is not a simple smooth function. 

8 This behavior is well-known from the single-impurity Anderson model [30]. 
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is seen to be enhanced. In particular, the momentum distribution 



1 f 



o 



dco ImG fc (ft)) 



(48) 



— oo 



has a discontinuity at the Fermi surface, given by n k - — n k + = (m*/m) , where kp — 

F F 

k F ±0+. 

4. DYNAMICAL MEAN-FIELD THEORY FOR CORRELATED 

LATTICE FERMIONS 



Itinerant quantum mechanical models such as the Hubbard model and its generalizations 
are much more complicated than classical, Ising-type models. Generally there do not 
even exist semiclassical approximations for such models that might serve as a starting 
point for further investigations. Under such circumstances the construction of a mean- 
field theory with the comprehensive properties of the Weiss molecular field theory 
for the Ising model will necessarily be much more complicated, too. As discussed 
above there do exist well-known mean-field approximation schemes, e. g. Hartree-Fock, 
random-phase approximation, saddle-point evaluations of path integrals, decoupling of 
operators. However, these approximations do not provide mean-field theories in the spirit 
of statistical mechanics, since they are not able to provide a global description of a 
given model (e.g., the phase diagram, thermodynamics, etc.) in the entire range of input 
parameters. 

Here the limit of high spatial dimensions d or coordination number Z has again been 
extremely useful [15]. It provides the basis for the construction of a comprehensive 
mean-field theory for lattice fermions which is diagrammatically controlled and whose 
free energy has no unphysical singularities. The construction is based on the scaled 
Hamiltonian (26) and the simplifications in the many-body perturbation theory discussed 
in Sec. 3.2. There we saw that the local propagator G(co), i.e., the amplitude for an 
electron to return to a lattice site, and the local but dynamical self-energy £(&)) are the 
most important quantities in such a theory. Since the self-energy is a dynamical variable 
(in contrast to Hartree-Fock theory where it is merely a static potential) the resulting 
mean-field theory will also be dynamical and can thus describe genuine correlation 
effects such as the Mott-Hubbard metal-insulator transition. 

The self-consistency equations of this dynamical mean-field theory (DMFT) for cor- 
related lattice fermions can be derived in different ways. Nevertheless, all derivations 
make use of the fact that in the limit of high spatial dimensions Hubbard-type models 
reduce to a "dynamical single-site problem", where the J-dimensional lattice model is 
effectively described by the dynamics of the correlated fermions on a single site which 
is embedded in a "bath" provided by the other particles. In the following I will present 
two different derivations of the single-site action and the self-consistency equations of 
the DMFT. I start with the approach by Janis [31, 32] who generalized the coherent 
potential approximation (CPA) for disordered systems. Then I discuss today's standard 
derivation developed independently by Georges and Kotliar [33] which is based on the 
mapping of the lattice problem onto a self-consistent single-impurity Anderson model; 
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this approach was also employed by Jarrell [34]. The DMFT equations derived within 
the CPA approach and the single-impurity approach, respectively, are identical. Never- 
theless it is the Anderson-impurity formulation which was immediately adopted by the 
community since it makes contact with the theory of quantum impurities and Kondo 
problems; for a review see Ref. [35]. This is a well-understood branch of many-body 
physics for whose solution efficient numerical codes had been developed already in the 
1980's, in particular by making use of the quantum Monte-Carlo (QMC) method [36]. 

Although the single-impurity based derivation of the DMFT is now the standard 
method I will present both approaches since it is always instructive to derive a theory in 
more than one way. 

4.1. Construction of the DMFT by generalizing the coherent potential 

approximation 

The coherent potential approximation (CPA) is a well-known mean-field theory which 
was originally developed in the context of disordered systems [37-40]. To be specific 
let us consider Anderson's tight-binding Hamiltonian 9 with local ("diagonal") disorder 
[41] 

Ha = — % £ + ( 49 ) 

VZ {Ri,Rj) 

where V, is a random variable drawn from some distribution function P(Vj). The elec- 
trons described by (49) do not interact. Therefore we deal with the problem of a single 
particle moving through a random medium. (Since there is no spin dependence one can 
suppress the spin index altogether and simply work with a spinless fermion). The prob- 
lem is made complicated by the randomness. It requires one to calculate the average of 
a physical quantity X (which is a function of all site energies V ; ) with respect to P(Vi). 
In particular, the arithmetic average is obtained by 

(X)a V :=n fdV i P(V i )X(V h ...,V L ). (50) 
Ri J 

One may now proceed as follows: 

1. The actual random medium, given by the local potentials Vi, is thought to be 
replaced exactly by an (unknown) effective medium, described by a complex, 
frequency-dependent self-energy; this defines the self-energy. 

2. Since the effective medium is required to yield an exact description of the random 
medium, we may remove the medium at a site R , replace it by an actual potential V) 
and then demand that, upon averaging, the scattering caused by the perturbation of 
the medium due to Vi vanishes identically. The self-consistency condition expressed 
in the last step determines the previously unknown self-energy. 



y Here the kinetic energy is already scaled according to (19) to allow for the limit Z — > °° to be taken, in 
which case the CPA becomes exact (see Sec. 4.1.1.). 
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Let Gij(z) = G be the Green function of the electron in the random medium, with z as a 
complex frequency and Go as the unperturbed Green function; we suppress site-indices 
for the moment [42]. The Lippmann-Schwinger equation for G is given by 

G = G + G VG or G 1 G=l+VG. (51) 

We now introduce a self-energy ^ij(z) = E into (51), which plays the role of an unknown 
potential 

(G 1 -E)G=1 + (V-E)G, (52) 

and demand 

(G) av = (G 1 -E)- 1 . (53) 
Multiplication of (52) by (G) av yields 

G=(G) av +(G) av (V-E)G. (54) 

Here V — E is a new scattering potential, whose effect can be described by a T-matrix 
via 

G=(G) av +(G) av r(G) av , (55) 

where 

V-L 

T = ; tt^- (56) 

l-(V-E)(G) av 

Averaging of (55) yields 

(r> av = 0, (57) 

which is a self-consistent equation for E. Eq. (57) demands that E is determined in such 
a way that the scattering due to the perturbation V — E vanishes. If (57) could be solved 
exactly the entire problem would be solved. However, an exact solution is usually not 
possible (an exception is the Lloyd model where P(Vi) is given by a Lorentzian), so 
that an approximation has to be made to proceed further. At this stage CPA assumes the 
self-energy to be site-diagonal 

Z ij ((o) = Z(co)8 ij , (58) 

i. e., to be homogeneous. Eq. (58) is equivalent to a single-site approximation and 
corresponds to step 2 in the construction of the CPA described below (50). Since E(co) is 
homogeneous it is a k independent, but frequency dependent potential and thereby only 
adds to the frequency dependence of G , i. e., the averaged propagator is simply given 
by the unperturbed propagator with shifted frequency : 

(G) av = G«.(z-E( z )). (59) 

For all sites Rj the condition (57) therefore reduces to (7]) av = 0, i. e., 

Vi-L(z) 



= 0, (60) 
l-(V,-E( Z ))Gg( Z -E(z))/ w 
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FIGURE 9. (a) The random medium, described by local potentials V„ is replaced by an unknown, but 
exact, effective potential E. (b) By demanding the average scattering from a single site with potential 
Vj — E to vanish, E is determined self-consistently; from Ref. [11]. 



where 

— OO 

is the local propagator, with N(co) as the DOS of the unperturbed system. Eq. (60) 
implies that in the effective medium the average scattering by a single site ("impurity") 
vanishes. The single-site aspect underlying the CPA may therefore be visualized as 
shown in Fig. 9 [39]. 

The CPA and its results have many attractive features: 

(i) CPA is a non-perturbative, but very simple and self-consistent theory; 

(ii) it may be considered the best single-site approximation for the disorder problem as 
can be inferred from the above derivation; 

(iii) it has the so-called Herglotz-properties E(z) = L*(z*) and Im(z — < for z < 
[44], which implies that it has the correct analytic properties (positive DOS, etc); 

(iv) it leads to very good qualitative and even quantitative results for the one-particle 

properties of disordered systems. The latter is true even in dimensions d ~ 3 and for 
parameter values of the disorder strength and the impurity concentration where CPA 
cannot be linked to perturbation theory. These properties have made CPA the most 
widely used approximation scheme for disordered systems. 



4.1.1. CPA and the limit d — >■ °° 



The single-site aspect of the CPA outlined above and, in particular, the property (58), 
indicate that the CPA will become exact in the limit of high coordination number Z. 
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Indeed, by investigating the moments of the electronic DOS it was observed that 1/Z 
appears as a "hidden" small parameter, which governs the size of the corrections to 
the CPA moments [45]. However, by considering only the moments one cannot draw 
conclusions about the validity of the CPA itself. 10 Using the scaling of t in (19) one can 
show [46] that (58), and hence CPA itself, indeed becomes exact for Z — > °°, irrespective 
of the lattice structure. In other words, CPA solves the disorder Hamiltonian (49) exactly 
in the limit Z — > <*>. This finding explains why for Z < <*> the CPA can be so successful 
even for intermediate values of the disorder strength and impurity concentration, i. e., 
when perturbation theory in these parameters is no longer justified. In fact, we now see 
that there is an additional small parameter, namely 1/Z, which allows for a perturbation 
expansion that is independent of the values of the input parameters. In view of the 
existence of a small parameter 1/Z, CPA is seen to be a controlled mean-field theory. 
It is therefore not so surprising that the CPA often gives qualitatively and quantitatively 

correct results even in dimensions d ~ 3. 



4.1.2. Alternative derivation of CPA 

I will now show that the CPA can also be derived from a variational principle [43]. 
To this end the coherent potential, i. e., the self-energy, will be determined from the 
averaged free-energy functional Q. of the corresponding single-site problem [47]. This 
field-theoretical approach has the advantage that it can be generalized to treat interacting 
lattice models, such as Hubbard-type models and disorder models on the same basis 
(see Sec. 4.1.3). Furthermore, the physical idea behind this single-site theory is very 
transparent and may be explained in terms of the following simple picture [43]. To 
calculate the averaged free energy corresponding to a single site J?,- of the medium we 
have to determine the energy density 

W« = (62) 

where L is the number of lattice sites. To this end we consider the second step of the 
CPA strategy outlined below (50) and drawn schematically in Fig. 10: 

(i) we start from the homogeneous, effective medium with free energy density Q med /L; 

(ii) we remove the medium at site Ri, i. e., subtract a corresponding energy Q ; ; and 

(iii) replace it by a site with a bare potential VJ, i.e., add a corresponding averaged energy 
(^ are ) av ; so we have 

M^ = %^-^. + (a^e )av _ (63) 



10 In fact, in Ref. [45] the kinetic energy was not scaled according to t — > f*/VZ, but the bandwidth was 
kept constant, corresponding to classical scaling (t — » t*/Z). In this case the Z — » °° limit does not lead to 
CPA at all, but becomes trivial. 
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(«0 0>) (c) 

FIGURE 10. Steps to construct the exact averaged free energy functional in d = °°: (a) Homogeneous 
effective medium, (b) the medium is removed at site /?,, (c) the cavity in the medium at site /?,• is filled with 
the actual potential V,. The probability amplitude for a particle from the medium to be at site i?, (more 
precisely: the return amplitude) is given by the local propagator G„; from Ref. [11]. 



To calculate the contributions in (63) we make use of the trace-log formula for the 
free energy density which in the non-interacting case (superscript 0) reads [48] 



f = -Ttr^ln[G° k (ico n )]-\ 



(64) 



Here G° k is given by (32), co n = (2n + 1)kT are the Matsubara frequencies and the trace 
operation implies 



a n=—<x> 

where in the present problem the spin is unimportant. Using (59), i. e. 



we have 



med 



<G,(z)) av = G£(z-L(z)) 



= -rtr£ln \G° k (ico n -L(ico n ))] 



-l 



(65) 



(66) 



= -T tr J dEN(E) In [ico n - L(i(O n ) +fi-E] 



(67) 



and 



Q, = -rtrln 



G° u {i(On -E(ifife))l 1 = ~T tr In [(G u (i(O n ))^] \ 



(68) 



where G Q U is the local, unperturbed propagator needed to remove the medium at site /?;. 
Finally, (Q^ 6 ) av can be obtained from 



a bare = _ r ln jpbare 



(69) 
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where the local partition function is determined by the action 5| )are as [48] 

^bare = J DcDc * £xp j _ ^bare^ c * }] _ (70) 

Here c,c* are Grassmann (anti-commuting) variables [48, 49]. The action is given by 

= -trc:[(G i; (/(0„)) av ]- 1 c„ + trc:(^-E(/(0„))c„, (71) 

where c n = c{i(O n ), c* = c*(i(D n ), and the two terms correspond to the kinetic and 
potential energy, respectively. We note that the local potential at site /f ; is measured 
relative to the surrounding effective medium E(z'co, 3 ). Since (71) is the expression for a 
non-interacting system (bilinear dependence on c and c*) the functional intergral in (70) 
is trivial to perform and yields, together with the trace-log formula detA = exp(trlnA), 

lni2f are = trln{[(G (7 (/^)) av ]- 1 -^ + L(/^)} (72a) 

= trln^- 1 ^)-^), (72b) 
where we introduced an effective local propagator £f by 

y-\i(On) = [(G„(/^)) av ]- 1 +E(/(0„). (72c) 

The ^(ico n ) propagator describes the coupling between the medium and the site Rj. 
The averaged free energy, which is a functional of ^ _1 (E), then follows as 

(Q) av = -LTtr{J dEN(E) In [ ico n + E(z'ft>„) - E ] 

-ln(^- 1 (/w„)-E(/(0„)) + (ln(^- 1 (/(0„)-y i O) av }. (73) 
By taking the variational derivative of (73) w.r.t. E, i. e., using the stationarity condition 

m 

we obtain 



^ = (74) 



' 1 (75, 



y-\i<Q n )-Y.{i<Q n ) \y-\iG} n )-Vi 

Together with (66) and (72c) this equation is seen to be identical to the self-consi stent 
Eq. (60) for L(ioo n ). Given a value ^^ 1 (ico n ) we obtain E(z'(0„) from (75), which deter- 
mines a new value ^^ l (i(O n ) = [G®((i(D n ) — L(ioo n ))]~ l + L(ioo n ) and so on. Eq. (75) 
expresses particularly clearly the single-site aspect of the CPA, as well as the role of 
E(z'(0„) as a homogeneous effective potential that describes the effect of the original 
random medium in the averaged system. 
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4.1.3. Generalization of the CPA approach to interacting systems 



The CPA was extensively used in the 1970's to investigate disordered systems. It was 
also applied to interacting models, e. g. the Hubbard model, by first transforming the 
model (approximately) to a random alloy problem ("alloy analogy") [40]. 

A new approach to the CPA, which makes use of field-theoretical functional integral 
techniques in connection with explicit diagrammatic perturbation theory, was initiated 
by Janis [47, 50]. Thereby the range of applicability of the CPA was extended to inter- 
acting lattice systems (spin systems and itinerant systems). In this generalized single- 
site approach the free-energy functional can be derived in closed form. The derivation 
is based on the general scheme for the construction of conserving approximations by 
Baym [51], i. e., the free energy Q{E}, which is a functional of E, is written as 

/3Q{E} = 4>{G} -tr(EG) -tr ln[(G°) _1 — E]. (76) 

Here G^ ((o) is the free propagator, (32), and G is the full propagator which is deter- 
mined by 

The quantity 4>{G} is obtained from a self-consistent perturbation expansion in G for 
the self-energy [51], the latter being defined by 

SG= Z - (77b) 

Hence Q is a functional of E only. To define Q. unambiguously one has to impose 
the boundary condition 4> = for E = 0. The construction of Q.{L} amounts to the 
construction of the functional 

A{E,G} = ${G}-tr(EG). (78) 

In a single-site theory, where the self-energy is purely local, A is fully determined by 
E and the local part of G. In this case the construction of e~P A reduces to a single-site 
problem, which can be performed explicitly [47, 50]. 

Due to the insight gained from the investigation of the large-J limit for fermionic 
lattice systems one can now conclude that the generalized CPA approach [47, 50] 
becomes exact in d = °° [3 1] , just as the CPA for disordered systems and other single-site 
theories discussed so far become exact in this limit. In particular, this field-theoretical 
approach can be used to derive the exact free energy for fermionic models in d = °° 
[31, 32, 43]. This leads to a comprehensive, controlled mean-field theory even for 
interacting fermionic models which is conceptually identical to the CPA for disordered 
systems or to the Weiss theory for the Ising model. Of course, this theory is necessarily 
much more complicated in detail than the previous mean-field theories since we now 
deal with a dynamical single-site problem in a fermionic bath. 

The physical idea behind the approach is the same as that described in the last 
subsection in connection with the CPA. Let us consider the motion of a particle on 
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a lattice in d = °°. The interaction with the other particles affects the motion. This 
change is exactly described by a yet unknown complex, dynamical field Y, a (cd). Hence 
the original system with its bare interactions has been exactly replaced by an effective 
medium; the latter is simply a system of non-interacting, itinerant electrons moving in a 
complex, homogeneous coherent potential L a (co). 



4.1.4. Exact free energy functional for the Hubbard model in d = °° and 

the self-consistency equations 

We will now use the generalized CPA described above to construct an exact expression 
for the free energy of the Hubbard model in d = °° [3 1, 32, 43]. We proceed as in the case 
of disordered systems (see Sec. 4.1.2 and Fig. 10), with VJ- replaced by v ia = Uh~i- a , but 
we do not have to perform any impurity average now. The single- site free energy density 
Q./Lis given by 

^mpd _ .-. k...... (79) 



T 



Using the analog of (66) 



G lha {z) = Gl a {z-U{z)) 
the first two terms are given by (see (66), (68)) 



I med 



= -rtr£ln[G« a (/^-E a (/^))] 1 



(80) 



(81) 



a- = -T tr ln[G°- CT (/ft)„ -E a (/w„))]- 1 = —T tr ln[G i7i ff (ifife)] -1 . (82) 

The contribution Q}? MS , obtained by replacing the medium on site by the actual, bare 
interaction (here: the Hubbard interaction) is again given by (69) and (70), where the 
single-site action is now given by [31, 43] 



S 



bare 



- ti C* o jGii^(i(0 n )] C a .n 



+ 



U J JT4(T)c t (T)c|(T)c ; (T)-tr4 ; „E a (/(0„)c a! „ 



(83) 



This expression has the same form as (71), but the spin-dependence of the problem has 
now been taken into account explicitly (c n — > c a , n ) and the local one-particle potential 
Vi has been replaced by the Hubbard interaction between up and down- spins. 
In analogy with (72c) we now introduce an effective local propagator ^ a {i(o n ) by 



y-\i<Qn) = [G<j(i(Qn)]- 1 +L G (iCO n ), 



(84) 



where G a (i0) n ) = Gu )G {iG) n ) . The effective propagator ^(iOn) again describes the cou- 
pling between the medium and the interaction-site R{. Note, that sites do not communi- 
cate with one another but only via the effective medium. With (84) the action (83) takes 
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the form 



S^{c a ,c* a ;^- 1 } = -trc^ n tf a l (ico n )c a , n + U f dxc\{x)c y {x)c\{x)c^x). (85) 

J 

In (70) the partition function j^ bare is given by an integral over anticommuting 
Grassmann variables. It may be transformed into a conventional functional integral 
over real, commuting variables by rewriting the Hubbard interaction in (85) using the 
Hubbard-Stratonovich transformation 

exp j-£/ j dxc* t (x)c t (x)cl(x) Ci (x)\ = J D^expj-^ J dx \r) 2 {x) + ?{x) 



-/ V W{^(T)[4(T)c t (T)+c|(T)c ; (T)]-/r 7 (T)[4(T)c t (T)-c|(T)c;(T)]}j j.(86) 

This is equivalent to the standard operator identity 

= ^ [(ntt + mi) 2 ~ (% _ "4) 2 ] ( 8V ) 

for the Hubbard interaction, where the two terms on the right-hand side correspond 
to charge and spin fluctuations, respectively. In (86) the fluctuations are described by 
real fluctuating fields t, (t) and Tj (t), respectively. Now that the interaction problem has 
been rewritten in terms of non-interacting particles in the presence of infinitely many 
fluctuating fields, the integration over the Grassmann variables in the expression for the 
partition function can be performed explicitly, yielding 



where now 



oo 



V=— oo 



(88) 



(89) 



with (|) mn = £ m _„, (fj) mn = r\ m _ n and (% l ) mn = 8 mn [%(i(O n )] l . 

Making use of the relation (64) and combining the three contributions in (72), the 
total free energy is then found as 

Q = -L T tr | J d(0N((0) \n[i(O n + ^l-L -Q)]- H% 1 - L a ) j - L T In are . (90) 

While Q. was originally a functional of L a it is now understood as a functional of 
[^(Efj)] -1 . The stationarity condition (74) is then 



8Q. 



= 0, 



(91) 
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which leads to 



1 




[%(ico n )}- 1 -L o (i(0 n ) 



x expf-S^,^ 1 }] 



(92b) 




1 



) 



nnl T],t; 



(92c) 



with (X)^ = (fDr]D^Xexp[-S^ K ])/^ aK . The right-hand side of (92a) is nothing 
but the very definition of the local propagator G,-, i(7 (z(0„) = G a (i(O n ) in terms of the 
action Sj^ 6 . Eqs. (92), together with 



and the local action (85) provide an exact, self-consistent set of equations for L a (or Sf CT ) 
for the Hubbard model in d = °° [31-33]. In the paramagnetic phase two quantities are 
seen to play the most important role: the local propagator G G {(6) and the self-energy 
E(co); this is illustrated in Fig. 11. 

Although the interaction problem is mathematically much more complicated than 
the analogous expression (38) for the disorder problem without interactions it can, 
in principle, again be solved by iteration: for given we obtain L a from (92), 
which yields a new ^r 1 via (93), etc. The exact local propagator is then provided by 
Gn^a = (^o 1 ~ ^a) 1 • As in the disorder case G ;!;t7 is completely expressed in terms of 
effective, averaged quantities. 

The expression used in (92c) brings out particularly clearly the similarities and differ- 
ences between the (on-site) interaction problem and the analogous expression (75) for 
the disorder case without interactions: 

(i) on the left-hand side of (92) the self-energy E appears again as a (homogeneous) 
effective medium, which is obtained exactly from the original system by some averaging 
process; 

(ii) however, this average is very different in the two cases: in the disorder problem it 
involves an integration over the actual disorder potentials Vj with a given disorder dis- 
tribution P(Vi), while in (92) it demands an integration over infinitely many fluctuating 
(random) fields r\ , £ , which simulate the actual interaction; 

(iii) the latter integration leads to a highly non-trivial coupling of the energies, i. e., 
Matsubara frequencies (£> n (note, that this coupling even exists in the static limit, i. e., 
for rjy = <^v = for v ^ 0), while in the disorder problem the corresponding Eq. (75) 
is diagonal in the frequency. This shows clearly that, although the interaction between 
electrons on different lattice sites has been reduced to an interaction of electrons with 




(93) 
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Z=12 

FIGURE 11. In the limit Z — > °° the Hubbard model effectively reduces to a dynamical single-site 
problem, which may be viewed as a lattice site embedded in a dynamical mean field. Electrons may hop 
from the mean field onto this site and back, and interact on the site as in the original Hubbard model (see 
Fig. 1). The local propagator G(co) (i.e., the return amplitude) and the dynamical self-energy of the 
surrounding mean field play the main role in this limit. The quantum dynamics of the interacting electrons 
is still described exactly. 



a mean field, the dynamics of the latter interaction is still non-trivial. Once more we 
observe that the many-body nature of the Hubbard model survives even in d = °°, making 
an analytic evaluation of the local propagator G a from (92), (93) generally impossible. 

4.2. Construction of the DMFT as a self -consistent single-impurity 

Anderson model 

Following the presentation of Georges, Kotliar, Krauth and Rozenberg [52] the dy- 
namical mean-field equations will now be derived using the so-called cavity method. 
This derivation starts by removing one lattice site together with its bonds from the rest 
of the lattice. The remaining lattice, which now contains a cavity, is replaced by a par- 
ticle bath which plays the role of the dynamical mean field (see Fig. 12). So far the 
derivation and the underlying physical picture coincides with that of the CPA approach 
described in the previous section. Now comes a new, physically motivated idea: the bath 
is coupled, via a hybridization, to the cavity. The resulting problem then amounts to the 
solution of an effective single-impurity Anderson model where the degrees of freedom 
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Z(k, co) 




FIGURE 12. The DMFT can be viewed as the mapping of the cZ-dimensional Hubbard model with 
self-energy £(£, co) onto a single site described by a single-impurity Anderson model which hybridizes 
with the bath provided by the other electrons [52]. In this local problem the self-energy T.(co) is now k 
independent; after Ref. [53]. 



of the bath, represented by an appropriate hybridization function, have to be determined 
self-consistently. 

To be specific, let us consider the partition function in the grand canonical ensemble 



iT= [l\Dct Dc ia exp[-S{cl,c ia }}. 



(94) 



The action S{c* a , c ;(T } for the Hubbard model is given by 



rP 

S{c* ,Cia} = / dx 

Jo 



L 4t ( t ) ( ^ - M ) c ia (T) + £ tijC* ia (t) c ja (T) 



i<7 



+ £f/4(T)c it (T)4(T)c^(T) 



(95) 



where again as in Sec. 4.1.2 we use Grassman variables c* iG , c, a - We split the action S 
into three parts 

S = S + AS + S ( -°\ (96) 



where Sq is the part containing only variables on site 



S = dx 
Jo 



E4a(T)(^-M)cOcr(T)+t/^ t (T)cot(T)c^(T)co|(T) 



(97) 



AS contains the hoppings between site and other sites of the lattice i ^ 

rj8 



fP 

AS= / dxY [tioc* ia (x)coa(x)+t 0i cl a (x)c ia (x)) 



(98) 
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and the rest, denoted by S^°\ is the part of the action where the site and its bonds are 
removed, i.e., with i, j^O 



S(°)= f dx 
Jo 



£ C/a(T)(4l-Mha(T)+ £ %4x(t)<7«t(t) 



+ t/£4(T)c ;t (T)4(T)c ; |(T) 
#0 



(99) 



Let us now rewrite the partition function 3f as 

^ = Yl Dc Oo Dc Oo^w[-So{co a ,co a }} 

x / II Dc* a Dc ia exp[-S {0 \c^c ia }]exp[-AS{c* 0a ,c 0a ,c^c ia }] (100) 

and use the following ensemble average 

<*>(o) = ^o) / .11 Oc^toXexph^^.Cto}] (101) 

taken with respect to S^ ) (the action where the site i = is excluded), with 3?W being 
the corresponding partition function. Then the partition function reads 

& = 2>{0) f YlDc^Dcoa&xpl-Soic^coa}] 
J a 

x (exp[-AS{c* 0a ,co a ,cl a ,c ia }}) m . (102) 

In the next step we expand the second exponent with respect to the action AS. As a result 
we obtain a formally infinite series with all possible many-particle correlation functions, 
i.e., 

% = ^(0) J YlD4 a Dc 0o exp[-S {4 a ,c 0o }} (l - j* dx(AS(x)) {0) 

+ dTl Jo rfT 2(A5(T 1 )A5(T 2 )) (0) + (103) 

where we used AS = Jq dxAS(r). In the fermionic case only the correlation functions 
with equal number of c and c* are non-zero. The lowest order term is second order and 
reads 



dT2(AS'(Ti)AS'(T2))(0) 

1 ^ 



/ dx x [ Jt 2 (A5(ti)A5(t 2 )) 
jo Jo 

1 /"P /"P 

= — / dTi dX2^ £ [^OA:(cy c7 (Tl)c fec7 (T2)>(0)CO(j(Tl)co (7 (T2) 



+ tojtko(Cja(Xl)cl a (X2)) { o ) C* Qa (Xi)co a (X2) ]. (104) 
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Higher-order terms are obtained similarly. The above expression can be rewritten with 
the use of one-particle correlation function 

G%{xx - t 2 ) = -{T,c ja {x x )ct a {x 2 )) {Q) (105) 

and takes the form 

1 fP fP 

-j o dx x ^ Jt 2 (A5(t 1 )A5(t 2 )) (0) 
= -f dx x f Jt 2 £ £ tjotkoGfjx, - x 2 )cl a {x l )co a {x 2 ). (106) 

Higher-order terms can be written in a similar way with the use of n-particle correlation 
functions. 

A non-trivial limit d — > °° is obtained by scaling the hopping amplitudes tjj as de- 
scribed in Sec. 3.2. For example, in the second-order contribution to the partition func- 
tion, (106), the hopping amplitudes must be scaled with Z^ Rq ~ r ^/ 2 because the one- 
particle correlation functions are proportional to 1 /zll^ 0_R >H/ 2 as discussed in Sec 3.2. 
In the calculation of higher-order terms we find that all connected higher-order terms 

vanish at least as 0(\/Z). Consequently, in the Z — > °° limit only the contribution G^ , 

or disconnected contributions made of products of G^ a 's remain. Applying the linked- 
cluster theorem and collecting only connected contributions in the exponential function 
one obtains the local action 



Sloe — 



fP d fP 

/ dxY,4 a {x){ 1 --ii)coo{x) + U dxc* ot (x)c ot (xy oi (x)c oi (x) 

JO (j OT JO 

+ f dx x f 3 Jt 2 £ £ t* t* k0 Gf^(x 1 -x 2 )c* 0a (x l )c 0o (x 2 ) 

J0 J0 cr j^O 



, (107) 



where the rescaled hoppings are denoted with a star. Introducing the hybridization 
function 

Aa(Ti - x 2 ) = - £ ty%Gf a {x x - x 2 ), (108) 

and employing the free ("Weiss") mean-field propagator one can express the DMFT 
local action in the following form (here the site index i = is omitted for simplicity) 

rP rP 

S\ oc = - dxi dx 2 Yc* a (xi)^ a l (xi-x 2 )c a (x 2 ) 

J0 J0 a 

+ U [ P dxc* t (x)c t (x)cl(x) Ci (x), (109) 

J 

% 1 (ti - x 2 ) = - - 8 Tl ? 2 - A a (xi - x 2 ). 



where 



(110) 



373 



Finally, we need the relation between the Green function G^ a (x — x') where the site 
i = is removed and the full lattice Green function, i.e., 

G fjl = G ij° ~ G iO<y G 00o G 0j° > ( 1 1 1 ) 

which holds for a general lattice. 

In order to obtain the full solution of the lattice problem it is convenient to express 
the relation between the local Green function Goocr = G a and the dynamical ("Weiss") 
mean field <S~ 1 in the form of a Dyson equation 11 

[G (i(d n )r l =&-\iCD n )-T, a (iCD n ) = iCD n + ii-Ao(iCD n )-I,o(iCD n ). (112) 
Then the lattice Green function (in &-space) G k< j{i(On) is given by 

G ka (iC0n) = ~. r- (H3) 

After performing the so-called lattice Hilbert transform we recover the local Green 
function 

G„(,*) = lG 4o (,*0=I— -^^-y dl4) 
After analytic continuation to real frequencies the local ("& averaged") propagator reads 

G ka ((0) = 1 - . r . (115) 

It is very important to realize that although the DMFT corresponds to an effectively 
local problem, the propagator G k {(o) is a momentum-dependent quantity. Namely, it 
depends on the momentum through the dispersion e k of the non-interacting electrons, 
but there is no additional momentum-dependence through the self-energy, since this 
quantity is strictly local within the DMFT. 

The set of self-consistent equations (109), (1 10), (112), (114) can be solved iteratively. 
In each step one solves the single-impurity problem given by the action (109), then one 
finds the new self-energy from the Dyson equation (112) and the new dynamical mean 
field from (114) and (112). The single-impurity problem is still a complicated many- 
body interacting problem which cannot, in general, be solved exactly. 



4.3. Solution of the DMFT self-consistency equations 

The dynamics of the full Hubbard model, (4), was found to remain complicated even 
in the limit d — > °° because of the purely local nature of the interaction. Hence an 



11 It should be noted that, in principle, any one of the local functions T. a (ico n ), Sf^ 1 (*'<»>,), or A CT (j'(8y,) can 
be viewed as a "dynamical mean field" acting on particles on a site, since they all appear in the bilinear 
term of the local action (109). 
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exact, analytic evaluation of the self-consistent set of equations (92), (93) for the local 
propagator G a or the effective propagator ^ a (i(D n ) is not possible. Exact evaluations are 
only feasible when there is no coupling between the frequencies. This is the case, for 
example, in the Falicov-Kimball model [54], which was solved analytically by Brandt 
and Mielsch [21] soon after the introduction of the d — > °° limit [15]. A valuable semi- 
analytic approximation is provided by the so-called iterated perturbation theory (IPT) 
[33], [55], [52]. 

Solutions of the general DMFT self-consistency equations require extensive numer- 
ical methods, in particular quantum Monte Carlo techniques [34], [56, 57], [52], the 
numerical renormalization group [58, 59], exact diagonalization [60-62], [52] and other 
techniques, whose discussion requires a separate series of lectures; here I refer the reader 
to the reviews quoted above. 

It quickly turned out that the DMFT is a powerful tool for the investigation of elec- 
tronic systems with strong correlations. It provides a non-perturbative and thermody- 
namically consistent approximation scheme for finite-dimensional systems which is par- 
ticularly valuable for the study of intermediate-coupling problems where perturbative 
techniques fail [63], [52, 64], [53], [65, 66]. 

In the remaining part of these lecture notes I shall discuss several applications of the 
DMFT to problems involving electronic correlations. In particular, I will address the 
Mott-Hubbard metal-insulator transition, and explain the connection of the DMFT with 
band-structure methods — the LDA+DMFT scheme — which is the first comprehensive 
framework for the ab initio investigation of correlated electron materials. 

5. THE MOTT-HUBBARD METAL-INSULATOR TRANSITION 

The correlation induced transition between a paramagnetic metal and a paramagnetic 
insulator, referred to as "Mott-Hubbard metal-insulator transition (MIT)", is one of the 
most intriguing phenomena in condensed matter physics [67-69]. This transition is a 
consequence of the competition between the kinetic energy of the electrons and their 
local interaction U . Namely, the kinetic energy prefers the electrons to move (a wave 
effect) which leads to doubly occupied sites and thereby to interactions between the 
electrons (a particle effect). For large values of U the doubly occupied sites become 
energetically very costly. The system may reduce its total energy by localizing the elec- 
trons. Hence the Mott transition is a localization-delocalization transition, demonstrating 
the particle-wave duality of electrons. 

Mott-Hubbard MITs are, for example, found in transition metal oxides with partially 
filled bands near the Fermi level. For such systems band theory typically predicts metal- 
lic behavior. The most famous example is V2O3 doped with Cr [70-72]. In particular, 
in (Vo.96Cro.o4)203 the metal-insulator transition is of first order below T = 380 K [71], 
with discontinuities in the lattice parameters and in the conductivity. However, the two 
phases remain isostructural. 

Making use of the half-filled, single-band Hubbard model (4) the Mott-Hubbard MIT 
was studied intensively in the past [5, 67-69, 72]. Important early results were obtained 
by Hubbard [73] within a Green function decoupling scheme, and by Brinkman and 
Rice [74] who employed the Gutzwiller variational method [4, 75], both at zero tem- 
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perature . Hubbard's approach yields a continuous splitting of the band into a lower 
and upper Hubbard band, but cannot describe quasiparticle features. By contrast, the 
Gutzwiller-Brinkman-Rice approach (for a review see Ref. [77]) gives a good descrip- 
tion of the low-energy, quasiparticle behavior, but cannot reproduce the upper and lower 
Hubbard bands. In the latter approach the MIT is signalled by the disappearance of the 
quasiparticle peak. 

To solve this problem the DMFT has been extremely valuable since it provided 
detailed insights into the nature of the Mott-Hubbard MIT for all values of the interaction 
U and temperature T [52, 78], [53]. 



5.1. DMFT and the three-peak structure of the spectral function 

The Mott-Hubbard MIT is monitored by the spectral function A((d) = — ^lmG(ft) + 
z'0 + ) of the correlated electrons 13 ; here we follow the discussion of Refs. [79], [53], [80]. 
The change of A(co) obtained within the DMFT for the one-band Hubbard model (4) at 
T = and half filling (n = 1) as a function of the Coulomb repulsion U (measured 
in units of the bandwidth W of non-interacting electrons) is shown in Figs. 13 and 14. 
While Fig. 13 is a schematic picture of the evolution of the spectrum when the interaction 
is increased, Fig. 14 shows actual numerical results obtained by the NRG [58, 80]. Here 
magnetic order is assumed to be suppressed ("frustrated"). 

While at small U the system can be described by coherent quasiparticles whose DOS 
still resembles that of the free electrons, the spectrum in the Mott insulator state consists 
of two separate incoherent "Hubbard bands" whose centers are separated approximately 
by the energy U. The latter originate from atomic-like excitations at the energies ±U/2 
broadened by the hopping of electrons away from the atom. At intermediate values of 
U the spectrum then has a characteristic three-peak structure as in the single-impurity 
Anderson model, which includes both the atomic features (i.e., Hubbard bands) and 
the narrow quasiparticle peak at low excitation energies, near co = 0. This corresponds 
to a strongly correlated metal. The structure of the spectrum (lower Hubbard band, 
quasiparticle peak, upper Hubbard band) is quite insensitive to the specific form of the 
DOS of the non-interacting electrons. 

The width of the quasiparticle peak vanishes for U — > U C 2(T). The "Luttinger pinning" 
at co = [20] is clearly observed. On decreasing U, the transition from the insulator to 
the metal occurs at a lower critical value U c \, where the gap vanishes. 

It is important to note that the three-peak spectrum originates from a lattice model 
with only one type of electrons. This is in contrast to the single-impurity Anderson 
model whose spectrum shows very similar features, but is due to two types of electrons, 



12 The Gutzwiller variational method [4, 75] consists of the choice of a simple projected variational wave 
function ("Gutzwiller wave function") and a semi-classical evaluation of expectation values in terms of 
this wave function ("Gutzwiller approximation"). As shown by Metzner and Vollhardt [76], [11] the 
Gutzwiller approximation becomes exact in the limit d — > °°. This result initiated the investigations of 
the Hubbard model in d — > °° [15] 

13 In the following we only consider the paramagnetic phase. 
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FIGURE 13. Evolution of the spectral function ("density of states") of the Hubbard model in the 
paramagnetic phase at half filling, a) non-interacting case, b) for weak interactions there is only little 
transfer of spectral weight away from the Fermi energy, c) for strong interactions a typical three-peak 
structure consisting of coherent quasiparticle excitations close to the Fermi energy and incoherent lower 
and upper Hubbard bands is clearly seen, d) above a critical interaction the quasiparticle peak vanishes 
and the system is insulating, with two well-separated Hubbard bands remaining; after Ref. [53]. 

namely the localized orbital at the impurity site and the free conduction band. Therefore 
the screening of the magnetic moment which gives rise to the Kondo effect in impurity 
systems has a different origin in lattice systems. Namely, as explained by the DMFT, 
the same electrons provide both the local moments and the electrons which screen these 
moments [52]. 

The evolution of the spectral function of the half-filled frustrated Hubbard model at 
finite temperatures, T = 0.0276 W, is shown in Fig. 15. This temperature is above the 
temperature of the critical point so that there is no real transition but only a crossover 
from a metallic-like to an insulating-like solution. The height of the quasiparticle peak 
at the Fermi energy is no longer fixed at its zero temperature value. This is due to a finite 
value of the imaginary part of the self-energy. The spectral weight of the quasiparticle 
peak is seen to be gradually redistributed and shifted to the upper (lower) edge of the 
lower (upper) Hubbard band. The inset of Fig. 15 shows the U -dependence of the value 
of the spectral function at zero frequency A ( £0 = 0). For higher values of U the spectral 
density at the Fermi level is still finite and vanishes only in the limit U — > °° (or for 
T -> 0, provided that U > U c2 (T = 0)). 

For the insulating phase DMFT predicts the filling of the Mott-Hubbard gap with 
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FIGURE 14. Evolution of the T = spectral function of the one-band Hubbard model with a semi- 
elliptic ("Bethe") DOS for interaction values U/W = 0,0.2,0.4, . . . , 1.6 (W: band width) calculated with 
the numerical renormalization group. At the critical interaction U&fW ~ 1.47 the metallic solution 
disappears and the Mott gap opens; from Ref. [80]. 



increasing temperature. This is due to the fact that the insulator and the metal are 
not distinct phases in the crossover regime, implying that the insulator has a finite 
spectral weight at the Fermi level. This behavior has been detected experimentally by 
photoemission experiments [81]. 

Altogether, the thermodynamic transition line U C (T) corresponding to the Mott- 
Hubbard MIT is found to be of first order at finite temperatures, being associated with a 
hysteresis region in the interaction range U c \ < U < U C 2 where U c \ and U C 2 are the val- 
ues at which the insulating and metallic solution, respectively, vanishes [52, 58], [82], 
[78, 79, 83], [53, 80]. The high precision, state-of-the-art MIT phase diagram by Bliimer 
[78] is shown in Fig. 16. The hysteresis region terminates at a critical point. For higher 
temperatures the transition changes into a smooth crossover from a bad metal to a bad 
insulator. 

It is interesting to note that the slope of the phase transition line is negative down to 
T = 0, which implies that for constant interaction U the metallic phase can be reached 
from the insulator by decreasing the temperature T, i.e., by cooling. This anomalous 
behavior (which corresponds to the Pomeranchuk effect [84] in 3 He, if we associate solid 
3 He with the insulator and liquid 3 He with the metal) can be easily understood from the 
Clausius-Clapeyron equation dU /dT = AS /AD. Here AS is the difference between the 
entropy in the metal and in the insulator, and AD is the difference between the number 
of doubly occupied sites in the two phases. Within the single-site DMFT there is no 
exchange coupling J between the spins of the electrons in the insulator, since the scaling 
(19) implies J oc —t 2 /U l/d — > for d — > °°. Hence the entropy of the macroscopically 
degenerate insulating state is 5i ns = k B In 2 per electron down to T = 0. This is larger than 
the entropy 5 me t °= T per electron in the Landau Fermi-liquid describing the metal, i.e., 
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FIGURE 15. Spectral function for the half-filled Hubbard model for various values of U at T = 
0.0276 W in the crossover region. The crossover from the metal to the insulator occurs via a gradual 
suppression of the quasiparticle peak at (0=0. The inset shows the U dependence of A(a>=0), in particular 
the rapid decrease for U w 1 . 1 W; from Ref. [79]. 



AS 

= ^met — >$ins < 0. At the same time the number of doubly occupied sites is lower 
in the insulator than in the metal, i.e., AD = D met — Di ns > 0. The Clausius-Clapeyron 
equation then implies that the phase-transition line T vs. U has a negative slope down to 
T = 0. However, this is an artifact of the single-site DMFT. Namely, there will always 
exist an exchange coupling between the electrons leading to a vanishing entropy of the 
insulator at T = 0. Since the entropy of the insulator vanishes faster than linearly with the 
temperature, the difference AS = S me t — Si ns eventually becomes positive, whereby the 
slope also becomes positive at lower temperatures 14 ; this is indeed observed in cluster 
DMFT calculations [85]. Since AS = at T = the phase boundary must terminate at 
T = with infinite slope. 

At half filling and for bipartite lattices in dimensions d > 2 (in d = 2 only at T = 0), the 
paramagnetic phase is unstable against antiferromagnetic long-range order. The metal- 
insulator transition is then completely hidden by the antiferromagnetic insulating phase, 
as shown in Fig. 17. 



Here we assume for simplicity that the metal remains a Fermi liquid, and the insulator stays param- 
agnetic, down to the lowest temperatures. In fact, a Cooper pair instability will eventually occur in the 
metal, and the insulator will become long-range ordered, too. In this case the slope dU/dT can change 
sign several times depending on the value of the entropy of the two phases across the phase transition. 
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FIGURE 16. High precision Mott-Hubbard MIT phase diagram showing the metallic phase and the 
insulating phase, respectively, at temperatures below the critical end point, as well as a coexistence region; 
fromRef. [78]. 
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FIGURE 17. On bipartite lattices and for half filling (n — 1) the paramagnetic phase is unstable against 
antiferromagnetism. The metal-insulator transition is then completely hidden by the antiferromagnetic 
insulating phase; from Ref. [86]. 

6. ELECTRONIC CORRELATIONS AND DISORDER 

The properties of real materials are strongly influenced by the electronic interaction and 
randomness. In particular, Coulomb correlations and disorder are both driving forces 
behind metal-insulator transitions (MITs) connected with the localization and dereal- 
ization of particles. While the Mott-Hubbard MIT is caused by the electronic repulsion 
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[5, 68, 87], the Anderson MIT is due to coherent backscattering of non-interacting parti- 
cles from randomly distributed impurities [41, 88]. The interplay between disorder and 
interaction effects gives rise to many fascinating phenomena [88-93], some of which 
will now be discussed following the presentation of Refs. [94-96]. 

The Mott-Hubbard MIT is characterized by the opening of a gap in the density of 
states at the Fermi level. By contrast, at the Anderson localization transition the character 
of the spectrum at the Fermi level changes from a continuous spectrum to a dense point 
spectrum. Both MITs can be characterized by a single quantity, the local density of states 
(LDOS). Although the LDOS is not an order parameter associated with a symmetry- 
breaking phase transition, it discriminates between a metal and an insulator, which is 
driven by correlations and disorder. 



6.1. Arithmetic vs. geometric averaging over the disorder 

The theoretical investigation of disordered systems requires the use of probability dis- 
tribution functions (PDFs) for the random quantities of interest. In physical or statistical 
problems one is usually interested in "typical" values of random quantities which are 
mathematically determined by the most probable value of the PDF, i.e., where the PDF 
becomes maximal. In many cases the complete PDF is not known, i.e., only limited 
information about the system provided by certain averages (moments or cumulants) is 
available. In this situation it is of great importance to choose the most informative av- 
erage of a random variable. For example, if the PDF of a random variable has a single 
peak and fast decaying tails this variable is usually well estimated by its first moment, 
known as the arithmetic average. The arithmetic average of a function F(e,-) is defined 
by (see also (50)) 



However, there are many examples, e.g., from astronomy, the physics of glasses or net- 
works, economy, sociology, biology or geology, where the knowledge of the arithmetic 
average is insufficient since the PDF is so broad that its characterization requires in- 
finitely many moments. Such systems are said to be non-self-averaging. One example is 
Anderson localization: when a disordered system is near the Anderson MIT [41], most 
of the electronic quantities fluctuate strongly and the corresponding PDFs possess long 
tails [97-100]. At the Anderson MIT the corresponding moments might not even exist. 
This is well illustrated by the local density of states (LDOS) of the system. The arith- 
metic mean of this random one-particle quantity does not resemble its typical value at 
all. In particular, it is non-critical at the Anderson transition [101] and hence cannot help 
to detect the localization transition. In this case the geometric mean [102-104], 



gives a much better approximation of the most probable ("typical") value of the LDOS. 
It vanishes at a critical strength of the disorder and hence provides an explicit criterion 
for Anderson localization [41, 105-107], [94-96]. 




(116) 



^geom = exp [(lnF(£f)) av ] , 



(117) 
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6.2. The Anderson-Hubbard model 



The fundamental electronic correlation model investigated here is the Anderson- 
Hubbard model 

H = -t J^c+c^ + ^eima + U^n^. (118) 

ia i 

The ionic energy e ; is a random, independent variable which describes the local, 
quenched disorder affecting the motion of the electrons. The disorder part is modeled 
by a corresponding PDF P(£i). For P(e,-) = the system is called pure. Here we use the 
continuous PDF 

P(ed= y \ — , (H9) 

with © as the step function. The parameter A is a measure of the disorder strength. 

A non-perturbative theoretical framework for the investigation of correlated lattice 
electrons with a local interaction is again given by the DMFT. If in this approach the 
effect of local disorder is taken into account through the arithmetic mean of the LDOS 
[108] one obtains, in the absence of interactions, the well-known coherent potential 
approximation [46] (see Sec. 4.1.), which does not describe the physics of Anderson 
localization. To overcome this deficiency Dobrosavljevic and collaborators formulated 
a variant of the DMFT where the geometrically averaged LDOS is computed from 
the solutions of the self-consistent stochastic DMFT equations [105] which is then 
incorporated into the self-consistency cycle [106]. Thereby a mean-field theory of 
Anderson localization can be derived which reproduces many of the expected features 
of the disorder-driven MIT for non-interacting electrons [106]. This scheme uses only 
one-particle quantities and is therefore easily incorporated into the DMFT for disordered 
electrons in the presence of phonons [109], or Coulomb correlations. In particular, the 
DMFT with geometrical averaging allows one to compute the phase diagram for the 
Anderson-Hubbard model with the continuous probability distribution function (1 19) at 
half filling [94]. In this way we find that, although the metallic phase is enhanced for 
small and intermediate values of the interaction and disorder, metallicity is eventually 
destroyed upon further increase of the disorder. Surprisingly, the Mott and Anderson 
insulators are found to be continuously connected. The phase diagram for the non- 
magnetic ground state is shown in Fig. 18. 

Interacting quantum many-particle systems with disorder pose fundamental chal- 
lenges for theory and experiment not only in condensed matter physics [68, 87-91], 
but most recently also in the field of cold atoms in optical lattices [110-115]. Indeed, 
ultracold gases have quickly developed into a fascinating new laboratory for quantum 
many-body physics [110, 111, 116-120]. A major advantage of cold atoms in optical 
lattices is the high degree of controllability of the interaction and the disorder strength. 
In particular, these quantum many-body systems will allow for the first experimental in- 
vestigation of the simultaneous presence of strong interactions and strong disorder. This 
very interesting parameter regime is not easily accessible in correlated electron mate- 
rials. Namely, at or close to half filling where interaction effects become particularly 
pronounced, strong disorder implies fluctuations (e.g., of local energies) of the order 
of the band width, which usually leads to structural instabilities. These limitations are 
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FIGURE 18. Non-magnetic ground-state phase diagram of the Anderson-Hubbard model at half filling 
as calculated by DMFT with the geometrically averaged ("typical") LDOS; after Ref. [94]. 



absent in the case of cold atoms in optical lattices where disorder can be tuned to be- 
come arbitrarily strong without destroying the experimental setup. Since at half filling 
and in the absence of frustration effects interacting fermions order antiferromagnetically, 
several basic questions arise: 

(i) How is a non-interacting, Anderson localized system at half filling affected by a local 
interaction between the particles? 

(ii) How does an antiferromagnetic insulator at half filling respond to disorder which in 
the absence of interactions would lead to an Anderson localized state? 

(iii) Do Slater and Heisenberg antiferromagnets behave differently in the presence of 
disorder? 

In Ref. [95] answers to the above questions were obtained by calculating the zero tem- 
perature, magnetic phase diagram of the disordered Hubbard model at half filling using 
DMFT with a geometric average over the disorder and allowing for a spin-dependence of 
the density of states (DOS). The results are collected in Fig. 19. Depending on whether 
the interaction U is weak or strong the response of the system to disorder is found to be 
very different. At strong interactions, U /W > 1, there exist only two phases, an AF insu- 
lating phase at weak disorder, A/W < 2.5, and a paramagnetic Anderson-Mott insulator 
at strong disorder, A/W > 2.5. The transition between these two phases is continuous. 
By contrast, the non-magnetic phase diagram for weak interactions, U /W < 1, has a 
much richer structure (Fig. 18). In particular, for weak disorder a paramagnetic metallic 
phase is found to be stable. It is separated from the AF insulating phase at large U by a 
narrow region of AF metallic phase. The AF metallic phase is long-range ordered, but 
there is no gap since the disorder leads to a redistribution of spectral weight [95]. 
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FIGURE 19. Magnetic ground-state phase diagram of the Anderson-Hubbard model at half filling as 
calculated by DMFT with the geometrically averaged ("typical") LDOS; after Ref. [95]. 



7. THEORY OF ELECTRONIC CORRELATIONS IN MATERIALS 

7.1. The LDA+DMFT approach 

Although the Hubbard model is able to explain basic features of the phase diagram of 
correlated electrons it cannot explain the physics of real materials in any detail. Clearly, 
realistic theories must take into account the explicit electronic and lattice structure of the 
systems. 

Until recently the electronic properties of solids were investigated by two essentially 
separate communities, one using model Hamiltonians in conjunction with many-body 
techniques, the other employing density functional theory (DFT) [121, 122]. DFT and 
its local density approximation (LDA) have the advantage of being ab initio approaches 
which do not require empirical parameters as input. Indeed, they are highly successful 
techniques for the calculation of the electronic structure of real materials [123]. How- 
ever, in practice DFT/LDA is seriously restricted in its ability to describe strongly cor- 
related materials where the on-site Coulomb interaction is comparable with the band 
width. Here, the model Hamiltonian approach is more general and powerful since there 
exist systematic theoretical techniques to investigate the many-electron problem with 
increasing accuracy. Nevertheless, the uncertainty in the choice of the model parameters 
and the technical complexity of the correlation problem itself prevent the model Hamil- 
tonian approach from being a flexible or reliable enough tool for studying real materials. 
The two approaches are therefore complementary. In view of the individual power of 
DFT/LDA and the model Hamiltonian approach, respectively, it had always been clear 
that a combination of these techniques would be highly desirable for ab initio investiga- 
tions of real materials, including, e.g., /-electron systems and Mott insulators. One of 
the first successful attempts in this direction was the LDA+U method [124, 125], which 
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combines LDA with a basically static, i.e., Hartree-Fock-like, mean-field approximation 
for a multi-band Anderson lattice model (with interacting and non-interacting orbitals). 
This method proved to be a very useful tool in the study of long-range ordered, insu- 
lating states of transition metals and rare-earth compounds. However, the paramagnetic 
metallic phase of correlated electron systems such as high-temperature superconductors 
and heavy-fermion systems clearly requires a treatment that goes beyond a static mean- 
field approximation and includes dynamical effects, e.g., the frequency dependence of 
the self-energy. 

Here the recently developed LDA+DMFT method, a new computational scheme 
which merges electronic band structure calculations and the dynamical mean-field the- 
ory, has proved to be a breakthrough [126-131], [53, 132-135]. Starting from conven- 
tional band structure calculations in the local density approximation (LDA) the corre- 
lations are taken into account by the Hubbard interaction and a Hund's rule coupling 
term. The resulting DMFT equations are solved numerically with a quantum Monte- 
Carlo (QMC) algorithm. By construction, LDA+DMFT includes the correct quasipar- 
ticle physics and the corresponding energetics. It also reproduces the LDA results in 
the limit of weak Coulomb interaction U. More importantly, LDA+DMFT correctly de- 
scribes the correlation induced dynamics near a Mott-Hubbard MIT and beyond. Thus, 
LDA+DMFT is able to account for the physics at all values of the Coulomb interaction 
and doping level. 

In the LDA+DMFT approach [126, 130, 131], [53] the LDA band structure is ex- 
pressed by a one-particle Hamiltonian H® DA , an( ^ * s men supplemented by the local 
Coulomb repulsion U and Hund's rule exchange J (here we follow the presentation of 
Ref. [80]). This leads to a material specific generalization of the one-band model Hamil- 
tonian 



H = Hl DA +U^n imt n imi + £ (V-8 aa ,J) 

n ima n im'a' ■ 

(120) 

m i i,m^m',a,a' 

Here m and m' enumerate the three interacting t2 g orbitals of the transition metal ion or 
the 4/ orbitals in the case of rare earth elements. The interaction parameters are related 
by V = U — 27 which holds exactly for degenerate orbitals and is a good approximation 
for the t2 g . The actual values for U and V can be obtained from an averaged Coulomb 
parameter U and Hund's exchange J, which can be calculated by constrained LDA. 

In the one-particle part of the Hamiltonian 



#lda = Hlda - ££ Ae rf w (121) 

m<7 



the energy term containing Ae^ is a shift of the one-particle potential of the interact- 
ing orbitals. It cancels the Coulomb contribution to the LDA results, and can also be 
calculated by constrained LDA [131]. 

Within the LDA+DMFT scheme the self-consistency condition connecting the self- 
energy E and the Green function G at frequency ft) reads: 

Gw m '(»)4 / d 3 k ([col + ill-H° DA (k)-IL(cQ)}- 1 ) (122) 

Vg J V / qm.q'm' 
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Here, 1 is the unit matrix, ji the chemical potential, H® DA (k) is the orbital matrix of 
the LDA Hamiltonian derived, for example, in a linearized muffin-tin orbital (LMTO) 
basis, E(ou) denotes the self-energy matrix which is nonzero only between the interacting 
orbitals, and implies the inversion of the matrix with elements n (=qm), n'{=q'm'), 
where q and m are the indices of the atom in the primitive cell and of the orbital, 
respectively. The integration extends over the Brillouin zone with volume Vg (we note 
that H® DA may include additional non-interacting orbitals). 

For cubic transition metal oxides Eq. (122) can be simplified to 



if the degenerate tig orbitals crossing the Fermi level are well separated from the other 
orbitals [131]. For non-cubic systems the degeneracy is lifted. In this case we employ 
Eq. (123) as an approximation, using different E m (co), N^(e) and G m {(6) for the three 
non-degenerate t% g orbitals. 

The Hamiltonian (120) is solved within the DMFT using standard quantum Monte- 
Carlo (QMC) techniques [36] to solve the self-consistency equations. From the imagi- 
nary time QMC Green function we calculate the physical (real frequency) spectral func- 
tion with the maximum entropy method [136]. 



7.2. Single-Particle Spectrum of Correlated Electrons in Materials 

Transition metal oxides are an ideal laboratory for the study of electronic correlations 
in solids. Among these materials, cubic perovskites have the simplest crystal structure 
and thus may be viewed as a starting point for understanding the electronic properties of 
more complex systems. Typically, the 3d states in those materials form comparatively 
narrow bands with width W ~ 2 — 3 eV, which leads to strong Coulomb correlations 
between the electrons. Particularly simple are transition metal oxides with a 3d 1 config- 
uration since, among others, they do not show a complicated multiplet structure. 

Photoemission spectra provide a direct experimental tool to study the electronic struc- 
ture and spectral properties of electronically correlated materials. Intensive experimental 
investigations of spectral and transport properties of strongly correlated 3d 1 transition 
metal oxides started with investigations by Fujimori et al. [137]. These authors observed 
a pronounced lower Hubbard band in the photoemission spectra (PES) which cannot be 
explained by conventional band structure theory. 

7.2.7. Excursion: Detection of electronic correlations in solids by 

photoemission spectroscopy 

In photoemission spectroscopy (PES) a photon of a given energy is used to emit 
an electron whose properties (energy, angular distribution) are measured in a detector. 
Angular resolved PES is referred to as ARPES. This technique measures occupied 
electronic states. This means that only the states described by the full spectral function 




(123) 
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FIGURE 20. Using photoemission spectroscopy (PES) the occupied electronic states can be measured. 
This corresponds to a multiplication of the (unknown) spectral function of a material with the Fermi 
function f((o) as shown in panel a). PES thus only measures the lower part of the full spectral function 
as shown in panel b). Inverse photoemission spectroscopy (IPES) or techniques like x-ray absorption 
spectroscopy (XAS) measure the unoccupied states. This corresponds to a multiplication of the spectral 
function with 1 — /(ft)) as shown in panel c). Hence IPES or XAS measures the upper part of the full 
spectral function as shown in panel d). 



of a material multiplied by the Fermi function /(ft)) are measured (see Fig. 20a), leading 
to the typical result shown in Fig. 20b. By contrast, inverse photoemission spectroscopy 
(IPES) measures the unoccupied electronic states. IPES is harder to perform and not as 
accurate as PES. But in many situations information about the unoccupied states is also 
available by X-ray absorption spectroscopy (XAS). In this case only the states described 
by the spectral function of a material multiplied by 1 — /(ft)) (Fig. 20c) are measured, 
leading to a typical result shown in Fig. 20d. 

Spectroscopic techniques are routinely used to investigate correlated electronic sys- 
tems. They provide very valuable information about the system since they measure the 
spectral function of a material which can be calculated theoretically. In particular, pho- 
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toemission techniques allow one to detect the correlation induced shift of spectral weight 
discussed in Sec. 5. 

In the following we will employ the LDA+DMFT approach to compute the k- 
integrated electronic spectra of two correlated materials, the 3d 1 system (Sr,Ca)V03 
and the charge-transfer insulator NiO. 



7.2.2. Sr x Cai- x V0 3 

SrV03 and CaVC>3 are simple transition metal compounds with a 3d 1 configuration 
(here we follow the presentation of Sekiyama et al. [138] and Nekrasov et al. [139]). 
The main effect of the substitution of Sr ions by the isovalent, but smaller, Ca ions is to 
decrease the V-O-V angle from 9 = 180° in SrV03 to 9 ps 162° in the orthorhombically 
distorted structure of CaV03. However, this rather strong bond bending results only in 
a 4% decrease of the one-particle bandwidth W and thus in a correspondingly small 
increase of the ratio U /W as one moves from SrV03 to CaV03. 

LDA+DMFT(QMC) spectra of SrV0 3 and CaV0 3 were calculated by Sekiyama et 
al. [138] by starting from the respective LDA DOS of the two materials; they are shown 
in Fig. 21. These spectra show genuine correlation effects, i.e., the formation of lower 
Hubbard bands at about 1.5 eV and upper Hubbard bands at about 2.5 eV, with well- 
pronounced quasiparticle peaks at the Fermi energy. Therefore both SrV03 and CaV03 
are strongly correlated metals. The DOS of the two systems shown in Fig. 21 are quite 
similar. In fact, SrV03 is slightly less correlated than CaV03, in accord with their 
different LDA bandwidths. The inset of Fig. 21 shows that the effect of temperature on 
the spectrum is small for T < 700 K. Spectra of SrV03 and CaV03 were also calculated 
independently by Pavarini et al. [140]. 

Since the three t2 g orbitals of this simple 3d 1 material are (almost) degenerate the 
spectral function has the same three-peak structure as that of the one-band Hubbard 
model shown in Fig. 15. The temperature induced decrease of the quasiparticle peak 
height is also clearly seen. As noted in Sec. 5 the actual form of the spectrum no 
longer resembles the input (LDA) DOS, i.e., it essentially depends only on the first 
three energy moments of the LDA DOS (electron density, average energy, band width). 
In the left panel of Fig. 22 the LDA+DMFT(QMC) spectra at 300K are compared with 
experimental high-resolution bulk PES. For this purpose the theoretical spectra were 
multiplied with the Fermi function at the experimental temperature (20 K) and Gauss 
broadened with the experimental resolution of 0.1 eV [138]. The quasiparticle peaks 
in theory and experiment are seen to be in very good agreement. In particular, their 
height and width are almost identical for both SrV03 and CaV03. The difference in 
the positions of the lower Hubbard bands may be partly due to (i) the subtraction of 
the (estimated) oxygen contribution which might also remove some 3d spectral weight 
below —2 eV, and (ii) uncertainties in the ab initio calculation of the local Coulomb 
interaction strength. In the right panel of Fig. 22 comparison is made with XAS data of 
Inoue et al. [141]. Core-hole life time effects were considered by Lorentz broadening the 
spectrum with 0.2 eV [142], multiplying with the inverse Fermi function (80K), and then 
Gauss broadening with the experimental resolution of 0.36eV [143]. Again, the overall 



388 




agreement of the weights and positions of the quasiparticle and upper ti g Hubbard band 
is good, including the tendencies when going from SrVC>3 to CaVC>3 (Cao.9Sro.1VO3 
in the experiment). For CaV03 the weight of the quasiparticle peak is somewhat lower 
than in the experiment. In contrast to one-band Hubbard model calculations, the material 
specific results reproduce the strong asymmetry around the Fermi energy w.r.t. weights 
and bandwidths. The results also give a different interpretation of the XAS than in Ref. 
[141] where the maximum at about 2.5 eV was attributed to an e g band and not to the ti g 
upper Hubbard band. The slight differences in the quasiparticle peaks (see Fig. 21) lead 
to different effective masses, namely m* /m = 2. 1 for SrV03 and m* jm = 2.4 for CaV03 . 
These theoretical values agree with m*/m = 2 — 3 for SrV03 and CaV03 as obtained 
from de Haas-van Alphen experiments and thermodynamics [144]. 

The experimentally determined spectra of SrV03 and CaV03 and the good agreement 
with parameter-free LDA+DMFT calculations confirm the existence of a pronounced 
three-peak structure in a correlated bulk material. Although the DMFT had predicted 
such a behavior for the Hubbard model (see Sec. 5.1.) it was not clear whether the 
DMFT result would really be able to describe real materials in three dimensions. Now 
it has been confirmed that the three-peak structure not only occurs in single-impurity 
Anderson models but also in three-dimensional correlated bulk matter. 



7.2.3. NiO 



Already in 1937, at the outset of modern solid state physics, de Boer and Verwey 
[145] drew attention to the surprising properties of materials with incompletely filled 
3J-bands, such as NiO. This observation prompted Mott and Peierls [146] to discuss 
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FIGURE 22. Comparison of the calculated, parameter-free LDA+DMFT(QMC) spectra of SrV0 3 
(solid line) and CaVC>3 (dashed line) with experiment. Left: Bulk-sensitive high-resolution PES (SrVC^: 
circles; CaV0 3 : rectangles). Right: Is XAS for SrV0 3 (diamonds) and Cao.9Sro.1VO3 (triangles) [141]. 
Horizontal line: experimental subtraction of the background intensity; after Ref. [139]. 



the interaction between the electrons. Ever since transition metal oxides (TMOs) were 
investigated intensively (here we follow the presentation of Runes et a/.[147]). This in- 
terest further increased when it was discovered that TMOs display an amazing multitude 
of ordering and electron correlation phenomena, including high temperature supercon- 
ductivity, colossal magnetoresistance and Mott metal-insulator transitions [1]. In the late 
1950's MnO and NiO were taken as the textbook examples of antiferromagnets. How- 
ever, when the importance of local Coulomb correlations in the transition metal J- shell 
was realized TMOs were considered candidates for Mott insulators [146]. In the mid 
1980's Zaanen, Sawatzky and Allen (ZSA) introduced their classification of TMOs and 
related compounds into Mott-Hubbard and charge-transfer (CT) systems [148]. In the 
early TMOs the ligand p-band is located well below the transition metal J-band and thus 
plays a minor role in the low energy dynamics. Such a case, called Mott-Hubbard sys- 
tem in the ZSA scheme, is well described by a multi-band Hubbard model. On the other 
hand, the late TMOs belong to the CT type where the p-band is situated between the 
interaction split d -bands. A more general Hamiltonian where the p- states are explicitly 
included is then needed, which can be viewed as a combination of multi-band Hubbard 
and Anderson lattice models. A major impulse for detailed investigations of CT systems, 
and especially of their hole doped regime, came with the discovery of high temperature 
superconductivity in cuprate perovskites. While the standard three-band Hamiltonian 
for cuprates [149] contains only one J-orbital per lattice site, the description of cubic 
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transition metal monoxides, the prominent member of which is NiO, requires the full 
set of d-orbitals. The latter are of interest not only for fundamental research, but play an 
important role also in fields such as geophysics [150]. Furthermore, recent progress in 
high pressure experiments [151] made the insulator-to-metal transition in some TMOs 
accessible in the laboratory, providing yet another stimulus for theoretical investigations. 

NiO is a type II antiferromagnet (7# = 523 K) with a magnetic moment of almost 
2\Xb and a large gap surviving well above T^. The standard LDA band theory predicts 
NiO to be a metal [152], or an antiferromagnetic insulator [153] if spin polarization 
is allowed. A severe underestimation of the gap and the magnetic moment suggests, 
however, that the Slater antiferromagnetic state obtained within LDA does not describe 
the true nature of NiO. On the other hand exact diagonalization studies on small clusters 
were quite successful in describing the single- and two-particle spectra [154], showing 
that the local Coulomb interactions are important. This made it clear that an explicit 
treatment of Coulomb interactions within the 3d shell is needed, and methods such as 
LDA+U [155], self-interaction correction [156], or GW [157] were applied. The static, 
orbitally dependent self-energy of LDA+U enforces a separation of the occupied and 
unoccupied d -bands and thus opens a gap comparable to experiment. This in turn leads to 
a significant improvement of the description of static properties such as the local moment 
or the lattice dynamics [158]. However, the LDA+U method is limited to an ordered 
state and does not yield the electronic excitations and the effect of doping correctly. A 
systematic inclusion of dynamical correlations is only made possible by the DMFT 

In Fig. 23 the calculated spectral densities [147] (resolved into Ni 3d and O 2p con- 
tributions) are compared to photoemission and inverse photoemission data [159]. Using 
the full p — d Hamiltonian it is possible to cover the entire valence and conduction bands 
spectra. Features corresponding to 45 and 4p bands at 10 eV and 13 eV, respectively, are 
not included in the theoretical spectrum. The relative intensity of the 2p contribution 
increases with decreasing photon energy [160]. Therefore the 120 eV spectrum is dom- 
inated by Ni 3d emission, while at 66 eV photon energy the O 2p contribution peaked 
around -4 eV is resolved (for a detailed orbital decomposition see Ref. [160]). The theo- 
retical spectrum very well reproduces the experimental features, including the size of the 
gap, the d character of the conduction band, the broad d peak at -9 eV, the position of the 
p-band, and the strong d contribution at the top of the valence band. While the gap and 
the Hubbard subbands can be described already with the static theory (LDA+U) [155], 
a dynamical treatment is apparently needed to capture the substantial redistribution of 
spectral weight between the incoherent (-9 eV) and resonant (-2 eV) features in the d 
spectrum. 

We see that by including the ligand p states and the on-site Coulomb interaction 
within the same framework one is able to provide a full description of the valence band 
spectrum and, in particular, of the distribution of spectral weight between the lower 
Hubbard band and the resonant peak at the top of the valence band. Good agreement 
with the available photoemission and inverse-photoemission data is found without need 
for adjustable parameters [147]. In a similar way ARPES data can also be explained in 
detail [161]. 
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FIGURE 23. Theoretical Ni-d (solid line) and O-p (shaded) resolved spectral densities compared to 
photoemission and inverse photoemission data obtained at 120 eV and 66 eV photon frequencies after 
Ref. [159]. Gaussian broadening of 0.6 eV full width at half maximum corresponding to the experimental 
resolution was applied to the theoretical curves. The circle marks position of the d 10 L excitation; from 
Ref. [147]. 



7.3. Correlation induced structural transformations 

In materials with correlated electrons the interaction between spin, charge, orbital, and 
lattice degrees of freedom leads to a wealth of ordering phenomena and complex phases 
[1]. The diverse properties of such systems and their great sensitivity with respect to 
changes of external parameters such as temperature, pressure, magnetic field or doping 
also make them highly attractive for technological applications [1]. In particular, orbital 
degeneracy is an important and often inevitable cause for this complexity [162]. A fasci- 
nating example is the cooperative Jahn-Teller (JT) effect — the spontaneous lifting of the 
degeneracy of an orbital state — leading to an occupation of particular orbitals ("orbital 
ordering") and, simultaneously, to a structural relaxation with symmetry reduction. 

Applications of LDA+DMFT so far mainly employed linearized and higher order 
muffin-tin orbital (L(N)MTO) methods [163] and concentrated on the study of correla- 
tion effects within the electronic system for a given ionic lattice (here we follow the pre- 
sentation of Leonov et al. [164, 165]). On the other hand, the interaction of the electrons 
with the ions also affects the lattice structure. LDA+DMFT investigations of particu- 
larly drastic examples, e.g., the volume collapse in paramagnetic Ce [166, 167] and Pu 
[168-170], and the magnetic moment collapse in MnO [171], incorporated the lattice by 
calculating the total energy of the correlated material as a function of the atomic volume. 
However, for investigations going beyond equilibrium-volume calculations, e.g., of the 
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FIGURE 24. Comparison of the total energies of paramagnetic KC11F3 computed by GGA and 
GGA+DMFT(QMC) as a function of the JT distortion. Error bars indicate the statistical error of the 
DMFT(QMC) calculations; from Ref. [164]. 



cooperative JT effect and other subtle structural relaxation effects, the L(N)MTO method 
is not suitable since it cannot determine atomic displacements reliably. This is partly due 
to the fact that the atomic-sphere approximation used in the L(N)MTO scheme, with 
a spherical potential inside the atomic sphere, completely neglects multipole contribu- 
tions to the electrostatic energy originating from the distorted charge density distribution 
around the atoms. 

Recently Leonov et al. [164, 165] formulated a computational scheme which now 
allows one to calculate lattice-relaxation effects caused by electronic correlations. To 
this end the GGA+DMFT — a merger of the "generalized gradient approximation" 
(GGA) and DMFT — was formulated within a plane-wave pseudopotential approach 
[172, 173]. Thereby the limitations of the L(N)MTO scheme in the direct calculation of 
total energies are overcome. In particular, this new method can be applied to determine 
the orbital order and the cooperative JT distortion in the paramagnetic phase of the 
prototypical JT system KCUF3. 

KCUF3 is long known to be a prototypical material with a cooperative JT distortion 
[162] where the electronic degrees of freedom are the driving force behind the orbital 
order [162, 174, 175]. Indeed, the relatively high (tetragonal) symmetry makes KCUF3 
one of the simplest systems to study. In particular, only a single internal structure 
parameter, the shift of the in-plane fluorine atom from the Cu-Cu bond center, is needed 
to describe the lattice distortion. 

The total energies as a function of the JT distortion 8jt = \ {di~ d s ) / (di + d s ) , where 
di and d s denote the long and short Cu-F bond distances, obtained by the GGA and 
GGA+DMFT, respectively, are compared in Fig. 24. It should be noted that the GGA 
not only predicts a metallic solution, but its total energy is seen to be almost constant 
for < 5jt < 4%. Both features are in contradiction to experiment since the extremely 
shallow minimum at Sjt — 2.5% would imply that KCUF3 has no JT distortion for T > 
100 K. By contrast, the inclusion of the electronic correlations among the partially filled 
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Cu e s states in the GGA+DMFT approach leads to a very substantial lowering of the total 
energy by ~ 175 meV per formula unit (fu). This implies that the strong JT distortion 
persists up to the melting temperature (> 1000 K), in agreement with experiment. The 
minimum of the GGA+DMFT total energy is located at the value 8jj = 4.2% which 
is also in excellent agreement with the experimental value of 4.4% [176]. This clearly 
shows that the JT distortion in paramagnetic KCUF3 is caused by electronic correlations. 

The GGA+DMFT scheme introduced in Refs. [164, 165] opens the way for fully 
microscopic investigations of the structural properties of strongly correlated electron 
materials such as lattice instabilities observed at correlation induced metal-insulator 
transitions. 

8. KINKS IN THE DISPERSION OF STRONGLY CORRELATED 

ELECTRON SYSTEMS 

The dispersion relation indicates at which energy and crystal momentum one-particle 
excitations can occur in a solid. The coupling between the excitations may lead to abrupt 
changes in the slope of the dispersion, referred to as "kinks". Such kinks thus carry 
important information about interactions in a many-body system (here we follow the 
presentation of Byczuk et al. [188, 189]). 

In systems with a strong electron-phonon coupling kinks in the electronic dispersion 
at 40-60 meV below the Fermi level are well known ("Kohn anomaly"). Therefore 
the kinks detected at 40-70 meV below the Fermi level in the electronic dispersion of 
high-temperature superconductors are taken as evidence for phonon [177, 178] or spin- 
fluctuation based [179, 180] pairing mechanisms. Kinks in the electronic dispersion at 
binding energies ranging from 30 to 800 meV are also found in various other metals 
[181, 182, 185-187] posing questions about their origins. 

Starting from the unexpected finding of kinks in the momentum-resolved spectral 
functions of SrV0 3 calculated by LDA+DMFT [190], Byczuk et al. [188] recently dis- 
covered a novel, purely electronic mechanism yielding kinks in the electron dispersion. 
This mechanism does not require a coupling of two different excitations as in previously 
known cases. The theory applies to strongly correlated metals whose spectral function 
shows well separated Hubbard subbands and central peak as, for example, in transition 
metal-oxides. For a microscopic description of these electronic kinks the single-band 
Hubbard model (4) with particle-hole symmetry was investigated by the DMFT at T = 
[188] as will be described below. 

The effective dispersion relation of the one-particle excitation is determined by 
the singularities of G(&, ft)) = (ft) + /I — £k — £(&, ft))) -1 , which give rise to peaks in the 
spectral function A ft)) = — lmG(k, Co)/k. Here ft) is the frequency, ji the chemical 
potential, the bare dispersion relation, and L(k, ft)) is the self-energy. If the damping 
given by the imaginary part of E(&, ft)) is not too large, the effective dispersion E k is thus 
determined by 

E k + n-£ k -ReZ(k,E k ) = 0. (124) 

Any kinks in E^ that do not originate from must therefore be due to changes in the 
slope of ReE(&, ft)) . 
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FIGURE 25. Local propagator and self-energy for a strongly correlated system, (a) Correlation-induced 
three-peak spectral function A(o) = —\mG{(o)/n with dips at ±£l = 0.45 eV. (b) Corresponding real 
part of the propagator, — ReG(ft)), with minimum and maximum at ±ft)max inside the central spectral 
peak, (c) Real part of the self-energy with kinks at ±0* (circles), located at the points of maximum 
curvature of ReG(ft)), (ft)* = 0.4ft> max = 0.03 eV). (d) ft) - 1/G(ft)) contributes to the self-energy. In 
general Re[co — 1/G(ft))] (full line) is linear in |ft)| < £1. The other contribution to the self-energy is 
— A(G(ft))) w — (m.2 — m\)G((o) (to lowest order in the moments m, of e^; here ni2 — m\=Q.5 eV 2 ). 
Therefore the nonlinearity of — Re[G(ft))] at ±0* determines the location of kinks; after Ref. [188]. 



The DMFT self-consistency equations are now used to express L(k,G)) = £(co) as 
L(co) = (0 + 11 — l/G((o) -A(G(c»)), where G((o) = J G(k,oo)dk is the fc-averaged 
local Green function and A(G) is an energy-dependent hybridization function, written 
here as a function of G{(o). 

Kinks in ReE(co) are found to appear at a new small energy scale which emerges 
quite generally for a three-peak spectral function A((o), see Fig. 25. Kramers-Kronig 
relations imply that Re[G(a))] is small near the dips of A{(o), located at +£l. Therefore 
Re[G(a>)] has a maximum and a minimum at iftWx inside the central spectral peak 
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(Fig. 25b). This directly leads to kinks in ReE(ft)) for the following reason. There are 
two contributions to E(co): 0) + jl - 1/G(ft)) and -A(G(ft))). While Re[ft) + /i - 1/G(ft>)] 
is linear in the large energy window \(o\ < £1 (Fig. 25d), the term — Re[A(G(ft)))] is 
approximately proportional to — Re[G(a>)] (at least to first order in a moment expansion), 
and thus remains linear only in a much narrower energy window \co\ < (Umax- The 
sum of these two contributions produces pronounced kinks in the real part of the self- 
energy at ±<o*, where ft)* = 0.41 (Om^x is the energy where Re[G(ft))] has maximum 
curvature (marked by circles in Fig. 25c). The Fermi-liquid (FL) regime with slope 
<9ReE(ft))/<9ft) = 1 — 1/Zfl thus extends only throughout a small part of the central 
peak (|ft)| < ft)*). At intermediate energies (ft)* < |ft)| < Q.) the slope is then given by 
<9ReE(ft))/<9ft) = 1 — 1/Zcp. The kinks at ±ft)* mark the crossover between these two 
slopes. As a consequence there is also a kink at ft)* in the effective band structure 
(Fig. 26). 

The FL regime terminates at the kink energy scale ft)*, which cannot be determined 
within FL theory itself. One of the most surprising results of the investigation is the fact 
that it is possible to describe properties of the system outside the Landau-Fermi-liquid 
regime fully analytically. Namely, one can express the quantities ft)* and Zqp in terms of 
Zfl and the bare density of states alone. Explicitly, one finds ft)* = OAlZp^D, where D 
is an energy scale of the noninteracting system, e.g., D is approximately given by half 
the bandwidth [188]. 

The energy scale ft)* involves only the bare band structure which can be obtained, for 
example, from band structure calculations, and the FL renormalization Zfl = 1/(1 — 
<9ReE(0)/<9ft)) = m/m* known from, e.g., specific heat measurements or many-body 
calculations. It should be noted that since phonons are not involved in this mechanism, 
ft)* shows no isotope effect. For strongly interacting systems, in particular close to a 
metal-insulator transition, ft)* can become quite small, e.g., smaller than the Debye 
energy. 

The theory described above explains the kinks in the slope of the dispersion as a direct 
consequence of the electronic interaction [188]. The same mechanism may also lead to 
kinks in the low-temperature electronic specific heat [191]. The kinks have also been 
linked to maxima in the spin susceptibility [192]. Of course, additional kinks in the 
electronic dispersion may also arise from the coupling of electrons to bosonic degrees 
of freedom, such as phonons or spin fluctuations. Interestingly, recent experiments [183] 
have found evidence for kinks in Ni(l 10), which maybe due to the electronic mechanism 
presented here. 

9. SUMMARY AND OUTLOOK 

Due to the intensive international research over the last two decades the DMFT has 
quickly developed into a powerful method for the investigation of electronic systems 
with strong correlations. It provides a comprehensive, non-perturbative and thermody- 
namically consistent approximation scheme for the investigation of finite-dimensional 
systems (in particular for dimension d = 3), and is particularly useful for the study of 
problems where perturbative approaches are inapplicable. For this reason the DMFT 
has now become the standard mean-field theory for fermionic correlation problems, in- 
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FIGURE 26. Kinks in the dispersion relation of the Hubbard model on a cubic lattice with interaction 
f/ =3.5 eV, bandwidth W ~ 3.46 eV, n= 1, implying a Fermi-liquid renormalization factor Zpp=0. 086. The 
intensity plot represents the spectral function A(k, co). Close to the Fermi energy the effective dispersion 
(white dots) follows the renormalized band structure = Zpp£jt (light line). For \ a>\ > (D* the dispersion 
has the same shape but with a different renormalization, E^ = Zcp£k — c sgn(£'^) (dark line). Here ©*=0.03 
eV, Zcp = 0.135, and c = 0.018 eV are all calculated from Zfl and (black line). A subinterval of T-R 
(white frame) is plotted on the right, showing kinks at ±o* (arrows); after Ref. [188]. 



eluding cold atoms in optical lattices [184], [120, 193]. The study of models in non- 
equilibrium within a suitable generalization of the DMFT has become yet another fasci- 
nating new research area [194-204]. 

Until a few years ago research into correlated electron systems concentrated on 
homogeneous bulk systems. DMFT investigations of systems with internal or external 
inhomogeneities such as thin films and multi-layered nanostructures are still very new 
[205-210]. They are particularly important in view of the novel types of functionalities 
of such systems, which may have important applications in electronic devices. Here the 
DMFT and its generalizations will certainly be very useful. 

In particular, the development of the ab initio band- structure calculation technique 
referred to as LDA+DMFT has proved to be a breakthrough in the investigation of 
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electronically correlated materials. It has already provided important insights into the 
spectral and magnetic properties of correlated electron materials, e.g., transition metals 
and their oxides [130, 131], [53, 132-134]. Clearly, this approach has a great potential 
for further developments. Indeed, it is not hard to foresee that the LDA+DMFT frame- 
work will eventually develop into a comprehensive ab initio approach which is able to 
describe, and even predict, the properties of complex correlated materials. 
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